NAVAL 

POSTGRADUATE 

SCHOOL 


MONTEREY,  CALIEORNIA 

THESIS 


CO-CHANNEL  INTERFERENCE  MITIGATION  USING 

SATELLITE  BASED  RECEIVERS 

by 

John  E.  Patterson 

December  2014 

Thesis  Advisor: 

Frank  E.  Kragh 

Second  Reader: 

Herschel  H.  Eoomis 

Second  Reader: 

Ric  Romero 

Approved  for  public  release;  distribution  is  unlimited 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


REPORT  DOCUMENTATION  PAGE 

Form  Approved  0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instruction, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington 
headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202—4302,  and 
to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  Blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

12-19-2014  Master’s  Thesis  07-02-2012  to  12-19-2014 

4.  TITLE  AND  SUBTITLE 

CO-CHANNEL  INTERFERENCE  MITIGATION  USING  SATELLITE  BASED  RE¬ 
CEIVERS 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

John  E.  Patterson 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

N/A 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

The  views  expressed  in  this  document  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the  Department  of 
Defense  or  the  U.S.  Government.  IRB  Protocol  Number:  N/A. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (maximum  200  words) 

This  thesis  investigates  a  means  to  mitigate  co-channel  interference  from  fourth  generation  cellular  signals  in  order  to  support  passive 
synthetic  aperture  radar  (SAR)  imaging  using  those  same  signals.  Synthetic  aperture  radar  is  a  staple  of  the  military  and  intelligence 
communities,  but  the  active  transmission  required  for  such  images  informs  the  target  of  the  imaging  process.  Use  of  passive  signals, 
such  as  the  orthogonal  frequency-division  multiplexing  (OFDM)  signals  of  fourth  generation  cellular  systems,  is  an  attractive  option, 
but  co-channel  interference  mitigation  is  required.  A  method  to  separate  the  transmitted  signals  that  leverages  the  estimated  signal 
delays  between  multiple  transmitters  and  receivers  is  examined  for  narrowband  and  wideband  signals.  Multiple  methods  of  assessing 
recovery  performance  are  proposed.  The  impact  of  noise  is  considered,  as  is  the  impact  of  collection  geometry  on  recovery  perfor¬ 
mance,  and  recovery  of  OFDM  signals  is  evaluated.  Signal  interpolation  is  a  critical  element  of  the  proposed  recovery  process,  and 
two  methods  are  compared  for  accuracy  and  speed  of  computation. 

14.  SUBJECT  TERMS 

Co-channel  interference,  orthogonal  frequency-division  multiplexing  (OFDM),  Passive  Synthetic  Aperture 
Radar,  Interpolation 

15.  NUMBER  OF 

PAGES  111 

16.  PRICE  CODE 

17.  SECURITY 

CLASSIFICATION  OF 

REPORT 

Unclassified 

18.  SECURITY 

CLASSIFICATION  OF  THIS 

PAGE 

Unclassified 

19.  SECURITY 

CLASSIFICATION  OF 

ABSTRACT 

Unclassified 

20.  LIMITATION  OF 
ABSTRACT 

UU 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  239-18 


1 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


CO-CHANNEL  INTERFERENCE  MITIGATION  USING  SATELLITE  BASED 

RECEIVERS 


John  E.  Patterson 
Commander,  United  States  Navy 
S.B.,  Massachusetts  Institute  of  Technology,  1998 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degrees  of 

ELECTRICAL  ENGINEER 

and 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  2014 


Author:  John  E.  Patterson 

Approved  by:  Prank  E.  Kragh 

Thesis  Advisor 

Herschel  H.  Eoomis 
Second  Reader 

Ric  Romero 
Second  Reader 

Clark  Robertson 

Chair,  Department  of  Electrical  and  Computer  Engineering 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


IV 


ABSTRACT 


This  thesis  investigates  a  means  to  mitigate  co-channel  interference  from  fourth  generation 
cellular  signals  in  order  to  support  passive  synthetic  aperture  radar  (SAR)  imaging  using 
those  same  signals.  Synthetic  aperture  radar  is  a  staple  of  the  military  and  intelligence 
communities,  but  the  active  transmission  required  for  such  images  informs  the  target  of  the 
imaging  process.  Use  of  passive  signals,  such  as  the  orthogonal  frequency-division  multi¬ 
plexing  (OFDM)  signals  of  fourth  generation  cellular  systems,  is  an  attractive  option,  but 
co-channel  interference  mitigation  is  required.  A  method  to  separate  the  transmitted  sig¬ 
nals  that  leverages  the  estimated  signal  delays  between  multiple  transmitters  and  receivers 
is  examined  for  narrowband  and  wideband  signals.  Multiple  methods  of  assessing  recovery 
performance  are  proposed.  The  impact  of  noise  is  considered,  as  is  the  impact  of  collection 
geometry  on  recovery  performance,  and  recovery  of  OFDM  signals  is  evaluated.  Signal 
interpolation  is  a  critical  element  of  the  proposed  recovery  process,  and  two  methods  are 
compared  for  accuracy  and  speed  of  computation. 
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Executive  Summary 


Introduction 

Imaging  of  the  earth  by  synthetic  aperture  radar  (SAR)  provides  a  number  of  advantages 
over  the  more  traditional  electro-optical  methods,  as  SAR  operates  without  a  requirement 
for  daylight  illumination  of  the  target  area  and  can  image  through  clouds,  providing  visi¬ 
bility  to  targets  that  would  not  otherwise  be  possible  to  image  [1].  Traditional  monostatic 
SAR  systems  require  illumination  of  the  target  area  with  a  transmitter,  which,  in  the  case 
of  military  or  intelligence  targets,  can  serve  to  alert  the  target  area  to  the  imaging  process. 
This  affords  the  opportunity  to  shelter  equipment  or  obscure  the  target  area  with  an  active 
radar  jammer,  negating  some  of  the  benefits  of  the  SAR  imaging  approach  [2]. 

Passive  radar  systems  exploit  signals  that  already  exist  in  the  target  area,  mitigating  the 
warning  provided  by  the  monostatic  SAR  system  [3].  The  proliferation  of  fourth  genera¬ 
tion  cellular  communications  networks  coupled  with  the  performance  of  their  orthogonal 
frequency-division  multiplexing  (OFDM)  waveform  makes  those  transmitters  an  appealing 
source  to  exploit  for  passive  SAR  imaging. 

Unfortunately,  that  same  proliferation  poses  difficulty  to  the  SAR  imaging  process  from 
the  satellite -based  receivers  in  the  form  of  co-channel  interference.  Cellular  networks  rely 
upon  spatial  frequency  reuse  in  order  to  make  the  best  possible  use  of  the  limited  spectrum 
available.  By  dividing  a  larger  geographic  area  into  smaller  cells  and  assigning  each  cell  a 
subset  of  the  available  frequencies,  it  is  possible  to  ensure  that  nearby  cells  do  not  reuse  a 
particular  frequency.  Satellite -based  receivers  view  many  of  these  transmitters  simultane¬ 
ously,  and  a  means  other  than  line-of-sight  is  required  to  isolate  the  individual  transmitter 
signals. 


Recovery  Process 

Consider  a  small  system,  consisting  of  two  transmitters  (transmitting  the  signals  5i(t)  and 
52 (t))  and  two  receivers  (which  receive  the  signals  ri  (t)  and  r2{t)).  In  that  case,  the  signals 


received  are  described  by 


nit)  =  Slit -Xn)+S2{t-Xi2)  (1) 

and 

r2it)  =  Siit-X2l)+S2it-X22),  (2) 

where  Xmn  is  the  delay  between  transmitter  n  and  receiver  m.  In  the  frequency  domain,  (1) 
and  (2)  become 

Riif)  =  (3) 

and 

Riif)  =  +S2(/)e-^'“"22,  (4) 

where  5„(/)  and  Rmif)  are  the  frequency-domain  representations  of  5„(t)  and  r^it),  re¬ 
spectively,  and  (O  =  2Ttf.  Generalizing  to  a  system  of  m  receivers  and  n  transmitters  gives 


.  .  Q  j^'^Xn 

'  Si(f)  ' 

_  RJf)  _ 

With  an  understanding  of  the  delays  x^n  between  each  transmitter  and  receiver,  it  is  possible 
to  solve  (5)  for  the  transmitted  signals  5i(/) . .  .5„(/).  It  is  important  to  recognize  that 
each  element  of  (5)  is  a  function  of  frequency  and  must  be  solved  for  each  frequency. 
Application  of  the  fast  Fourier  transform  (FFT)  to  the  sampled  time-domain  signals  yields 
signal  descriptions  in  the  discrete  frequency  domain  and  facilitates  this  solution. 


OFDM  Signal  Recovery 

A  basic  system  of  two  receivers  collecting  the  signals  of  two  OFDM  transmitters  is  ad¬ 
equate  to  demonstrate  the  potential  recovery  performance.  The  transmitter  and  receiver 
geometry  is  depicted  in  Figure  1,  and  the  corresponding  distances  used  to  compute  the  sig¬ 
nal  propagation  delays  are  provided  in  Table  1.  The  OFDM  signal  parameters  for  each  of 
the  transmitters  are  summarized  in  Table  2. 

xviii 


Table  1.  Transmitter  and  receiver  geometries.  (In  kilometers.) 


Element  x  y  alt 

Txl  0  0  (T 

Tx  2  0  30  0 

Rx  1  146.2  0  400 

Rx2  146.2  50  400 


Front  Top 


Figure  1.  Overview  of  test  geometry. 

Table  2.  OFDM  test  signal  parameter  summary. 

Value 

Subcarriers  6 

Modulation  Scheme  QPSK 
Bandwidth  (kHz)  600 

Symbols  16 


In  order  to  quantitatively  assess  recovery  performance,  the  correlation  coefficient 


Z(s,s) 


s*s 


(6) 


is  used  to  compare  the  original  signal  s  with  the  recovered  signal  s.  The  scalar  product 
of  two  signals  is  a  means  of  assessing  the  similarity  of  two  signals  [4].  The  correlation 
coefficient  normalizes  the  scalar  product  by  signal  powers,  such  that  x{)  varies  from  -1 
to  1,  where  a  value  of  1  indicates  proportional  signals,  0  indicates  orthogonality,  and  -1 
indicates  proportional  signals  with  different  sign  convention  [5]. 


The  correlation  coefficients  of  the  recovered  signals  are  provided  in  Table  3  and  indicate 
that  the  recovered  signals  are  very  closely  matched  to  the  original  signals.  A  graphical 


XIX 


comparison  of  a  subset  of  the  time-domain  samples  of  the  original  and  recovered  signals^ 
is  provided  in  Figure  2,  and  a  frequency-domain  comparison,  a  subset  of  the  total  signal 
spectrum,  is  provided  in  Figure  3.  Both  the  time  and  frequency  perspectives  indicate  excel¬ 
lent  recovery  performance.  The  frequency-domain  perspective  does  indicate  a  discrepancy 
at  the  /  =  0  point,  which  is  expected  as  the  delay  matrix  in  (5)  is  singular  at  that  point. 

Recovered  signal  error  is  computed  using 

Ss=S-S  (7) 

in  the  frequency  domain,  where  S  is  the  vector  of  original  signal  samples  and  S  is  a  vector 
of  the  recovered  signal  samples.  Similarly,  time-domain  error  is  computed  using 

ss  =  s-s,  (8) 

where  s  and  s  are  the  time-domain  samples  of  the  original  and  recovered  signals. 
Frequency-  and  time-domain  recovery  error  plots  are  provided  in  Figure  4.  Given  the 
known  singularity  at  /  =  0,  that  point  is  omitted  from  the  frequency-domain  error  plot. 
Both  error  plots  indicate  good  recovery  performance,  with  frequency-domain  errors  on  the 
order  of  computer  rounding  error.  The  steady-state  error  in  Figure  4b  is  likely  a  result  of 
the  singularity  at  /  =  0. 

Conclusion 

The  above  results  show  good  correlation  between  the  original  and  recovered  signals.  While 
additional  work  is  required  to  ensure  that  the  recovered  signals  sufficiently  isolate  desired 
signals  from  co-channel  interference,  the  method  appears  to  demonstrate  the  potential  to 
support  that  isolation  in  a  passive  SAR  system. 

Table  3.  Recovery  performance  figures. 

Correlation 

Signal  1  0.9995  - 4.6035  x 
Signal  2  0.9993  +  3.0037  x  lO-^^f 


'All  signals  are  simulated  and  assumed  to  be  either  voltage  or  current.  Units  are  unspecified  in  the  figure. 
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Figure  2.  Time  domain  of  transmitted  and  recovered  signals. 
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CHAPTER  1: 
Introduction 


Synthetic  aperture  radar  (SAR)  imaging  was  conceived  in  195 1  by  Carl  Wiley,  and  the  util¬ 
ity  of  the  technique  was  quickly  recognized.  The  early  work  in  the  field  depended  upon  the 
ability  to  finely  resolve  the  Doppler  characteristics,  or  changes  in  the  signal  return  due  to  the 
relative  motion  between  radar  and  target,  of  the  received  radar  signal  [1].  Initial  approaches 
relied  upon  optical  processing  techniques  and  utilized  a  variety  of  spherical,  conical  and 
cylindrical  lenses  to  create  the  radar  image  [1].  Advances  in  modern  computing  and  dig¬ 
ital  signal  processing  in  the  1980s  facilitated  the  more  modern  wavefront  reconstruction 
techniques  that  commonly  see  use  today  [2].  SAR  technological  developments  have  facil¬ 
itated  a  proliferation  of  SAR  imaging  systems  including  the  Shuttle  Imaging  Radar  (SIR- 
C/X-SAR)  and  Canada’s  RADARSAT  mission,  which  are  both  examples  of  SAR  satellites 
developed  for  Remote  Sensing  missions  studying  the  earth  [1].  Regardless  of  the  ultimate 
use,  SAR  provides  a  unique  means  to  monitor  the  earth  that  is  capable  of  overcoming  some 
considerable  limitations  of  the  more  commonly  known  electro-optical  systems.  As  an  ex¬ 
ample,  the  study  of  tropical  regions  by  remote  sensing  scientists  was  particularly  hampered 
by  cloud  cover  in  those  tropical  regions  as  even  after  a  decade  of  imaging  with  the  electro- 
optical  and  infrared  sensors  on  the  LANDSAT  satellites,  there  were  still  regions  of  the  earth 
that  had  not  been  imaged  without  cloud  cover  [1].  SAR  provides  a  means  to  collect  useful 
data  for  those  regions  despite  the  near  constant  cloud  coverage. 

1.1  Overview 

Much  like  the  benefits  afforded  to  the  remote  sensing  community,  SAR  offers  a  number  of 
advantages  over  electro-optical  imaging  systems  to  the  military  and  intelligence  commu¬ 
nities.  The  ability  to  image  a  target  area  when  obscured  by  clouds  is  of  equal  importance 
to  the  military  and  intelligence  communities.  The  ability  to  image  a  target  without  a  de¬ 
pendence  on  illumination  from  the  sun  allows  targets  to  be  monitored  at  all  hours  of  the 
day. 

Unfortunately,  and  unlike  the  imaging  provided  by  electro-optical  systems,  the  transmit¬ 
ter  required  for  traditional  monostatic  SAR  systems  reveals  the  imaging  process  when  that 
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signal  is  detected  in  the  target  area,  affording  the  target  the  opportunity  to  employ  coun¬ 
termeasures  to  prevent  the  imaging  of  sensitive  target  area  information.  In  addition  to 
sheltering  equipment,  the  use  of  an  active  jammer  when  the  SAR  transmitter  is  detected 
might  also  preclude  the  successful  imaging  of  the  target  area  [3]. 

By  exploiting  transmissions  that  already  exist  in  the  target  area,  a  passive  SAR  system  mit¬ 
igates  some  of  the  risks  associated  with  operations  targeting  an  adversary  aware  of  the  col¬ 
lection.  Gutierrez  del  Arroyo  [4]  highlights  a  number  of  potential  advantages  of  a  passive 
SAR  system,  including  both  covertness  and  decreased  susceptibility  to  electronic  attack. 
The  passive  exploitation  of  existing  signals  necessitates  a  radar  designed  to  operate  with 
its  receiver  well  separated  from  the  transmitter.  Radar  systems  that  exhibit  this  separation 
between  the  radar  transmitter  and  receiver  are  known  as  bistatic  radars  [5]. 

A  variety  of  signal  sources  might  be  considered  for  use  with  such  a  passive  SAR  sys¬ 
tem.  Two  commercially  developed  systems,  Lockheed  Martin’s  Silent  Sentry  and  Thales’ 
Homeland  Alerter  100,  utilize  analog  FM-radio  transmissions  for  the  detection  and  track¬ 
ing  of  aircraft.  These  are  the  first  two  known  commercial  systems  to  utilize  passive  radar 
techniques,  and  the  random  nature  of  the  FM  transmissions  that  they  utilize  as  illuminators 
provides  a  signal  processing  challenge  [4].  Orthogonal  frequency-division  multiplexing 
(OFDM)  is  a  digital  waveform  that  enjoys  wide  use  in  communications  applications.  It  has 
also  been  proposed  and  extensively  studied  as  a  potential  radar  waveform,  with  a  number 
of  characteristics  that  make  it  well  suited  to  the  radar  application  [4],  [6].  One  of  the  more 
prevalent  uses  of  the  OFDM  waveform  is  in  fourth  generation  cellular  communication  net¬ 
works  including  the  IEEE  802.16  Wireless  Metropolitan  Area  Networks  compliant  WiMax 
technology  and  long  term  evolution  (LTE)  [7],  [8].  The  proliferation  of  this  signal  coupled 
with  the  low  range  and  doppler  sidelobes  of  its  ambiguity  function  makes  it  an  excellent 
source  for  passive  SAR  imaging  [9],  [10]. 

That  is  not  to  say,  however,  that  using  OEDM-based  cellular  transmissions  is  without  chal¬ 
lenges.  Eundamental  to  cellular  technology  is  the  concept  of  spatial  frequency  reuse.  In 
order  to  support  a  large  number  of  users  with  a  limited  amount  of  available  spectrum,  the 
cellular  systems  must  make  use  of  the  same  frequencies  multiple  times  in  a  given  geo¬ 
graphic  area.  As  explained  in  [11],  a  given  coverage  area  is  subdivided  into  cells  (often 
modeled  as  hexagons,  but  much  more  irregular  in  reality).  See  Eigure  1.1  for  an  illustration 
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of  potential  cell  layouts.  Each  cell  is  assigned  a  subset  of  the  available  frequencies,  such 
that  no  adjacent  cells  utilize  the  same  frequencies.  A  frequency  reuse  distance  is  estab¬ 
lished  such  that  co-channel  interference  does  not  impact  the  communication  channel,  and 
only  outside  of  that  distance  from  a  cell  may  the  frequencies  be  used  again.  This  cellular 
allocation  of  frequencies  facilitates  an  increased  density  of  transmitters,  which  in  turn  al¬ 
lows  the  cell  phones  and  transmitters  to  operate  at  lower  power  levels.  This,  in  combination 
with  line-of-sight  limitations,  mitigates  interference  between  cells  [12]. 


Figure  1.1.  Frequency  reuse  in  adjacent  cells,  after  [11]. 

While  spatial  frequency  reuse  is  a  convenient  way  to  maximize  the  utilization  of  the  avail¬ 
able  spectrum,  it  poses  a  challenge  to  an  asset  attempting  to  use  such  a  signal  as  an  illu¬ 
minator  for  a  SAR  image.  The  SAR  image  reconstruction  process  utilizes  a  matched  filter 
and,  therefore,  requires  a  detailed  understanding  of  the  signal  used  as  the  illuminator  for 
the  image.  It  is  likely  that  the  signals  from  surrounding  transmitters  will  also  be  received, 
presenting  a  co-channel  interference  issue  that  must  be  mitigated  before  the  SAR  image 
can  be  formed.  That  mitigation  is  the  focus  of  this  effort. 

1.2  Literature  Review  and  Prior  Work 

Levanon  introduced  the  use  of  multi-carrier  phase  coded  signals  (now  more  commonly 
referred  to  as  OFDM)  and  presented  an  analysis  of  their  characteristics  as  a  radar  illumi¬ 
nator  [10].  Gutierez  del  Arroyo  and  Jackson  have  published  a  number  of  papers  related  to 
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the  use  of  eommereial  transmitters  as  illuminators  for  passive  SAR,  ineluding  the  devel¬ 
opment  of  a  bistatie  ambiguity  funetion  for  WiMax  emitters  [9],  eolleetion  and  proeessing 
of  WiMax  transmissions  for  bistatic  imaging  [6],  and  analysis  of  the  impact  of  WiMax- 
specific  signal  characteristics  [4].  That  work  assumes  that  interference  from  neighboring 
cells  is  mitigated  through  network  design,  which  is  reasonable  at  low  collection  altitudes 
but  becomes  problematic  with  collections  from  higher  altitude  given  the  increased  line-of- 
sight  to  numerous  cells  in  the  network. 

Blind  source  separation  (BSS)  is  a  field  dedicated  to  the  recovery  of  original  signals  when 
numerous  signals  have  been  combined  prior  to  receipt.  There  are  a  number  of  approaches 
to  signal  separation  that  might  be  applied  to  this  problem.  Stone  presents  an  overview 
of  one  BSS  approach  know  as  independent  component  analysis  (ICA)  [13].  While  some 
work  in  the  BSS  field  has  been  with  a  limited  number  of  receivers  [14],  the  bulk  of  that 
research  is  for  situations  with  equal  numbers  of  transmitters  and  receivers  as  well  as  for 
real-valued  signals,  often  from  audio  sources.  Douglas  has  extended  blind  separation  to 
complex  valued  signals  [15],  and  many  of  the  leading  algorithms  for  blind  source  separation 
are  collected  in  [16].  Those  techniques  often  depend  upon  the  non-Gaussianity  of  the 
original  signals,  and  OFDM  signals  were  not  readily  separated  in  initial  attempts  to  apply 
those  tools. 

Fason  [17]  and  Holder  [18]  developed  a  SAR  model  in  MATLAB  that  can  be  used  for 
analysis  of  various  parameters  that  impact  SAR  performance.  The  work  described  in  this 
thesis  is  intended  as  a  precursor  to  the  extension  of  that  model  to  the  bistatic  case  and 
performance  analysis  using  recovered  OFDM-based  signals. 

1.3  Outline 

An  overview  of  the  separation  approach  is  provided  in  Chapter  2,  including  a  description  of 
the  methods  used  for  calculating  the  necessary  signal  delays.  A  number  of  signal  recovery 
performance  parameters  are  proposed,  and  separation  is  evaluated  with  both  narrowband 
and  wideband  signals.  The  frequency  response  that  results  from  delayed  or  shifted  ver¬ 
sions  of  those  signals  is  evaluated.  Recovery  performance  under  noisy  conditions  is  also 
evaluated. 

A  description  of  the  OFDM  signal  model  used  is  provided  in  Chapter  3  along  with  an 
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analysis  of  the  recovery  of  a  number  of  OFDM  signals  with  various  constituent  parameters. 
That  work  is  extended  to  several  collection  geometries  in  Chapter  4,  and  a  means  to  evaluate 
potential  collection  geometries  prior  to  signal  collection  in  order  to  increase  the  likelihood 
of  viable  signal  recovery  is  also  discussed. 

Results  are  summarized  in  Chapter  5,  and  a  number  of  areas  requiring  further  study  are 
identified  and  discussed. 

An  analysis  of  two  different  methods  of  imposing  a  time  delay  on  a  sampled  signal  is  pro¬ 
vided  in  Appendix  A.  Figures  detailing  the  results  of  the  signal  recovery  efforts  in  Chapters 
3  and  4  are  provided  in  Appendix  B.  MATLAB  functions  that  implement  critical  elements 
of  the  work  described  throughout  this  thesis  are  provided  in  Appendix  C. 
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CHAPTER  2: 
Signal  Recovery 


Cellular  technologies  capitalize  on  line-of-sight  limitations  in  order  to  use  available  spec¬ 
trum  efficiently.  These  line-of-sight  limitations  are  negated  when  the  receiver  is  elevated, 
with  potential  to  receive  multiple  signals  simultaneously.  In  many  applications,  use  of  those 
signals  requires  the  ability  to  isolate  a  single  transmission  from  a  multitude  of  co-channel 
transmitters,  often  without  detailed  knowledge  of  the  transmitted  signals. 

The  main  contribution  of  the  work  described  in  this  thesis  is  a  method  to  separate  co¬ 
channel  signals  by  leveraging  the  known  positional  information  of  both  the  transmitters 
and  receivers  to  facilitate  the  recovery  of  the  desired  signal. 


2.1  General  Approach 

Given  a  transmitter  sending  a  signal  s{t),  a  receiver  will  receive  r{t)  =  s{t  —  t),  where 
T  is  the  delay  due  to  propagation  between  the  transmitter  and  receiver,  a  function  of  the 
distance  and  propagation  rate  of  the  intervening  medium,  and  attenuation  and  noise  effects 
are  neglected.  Considering  a  small  system  of  two  transmitters  and  two  receivers,  the  signals 
received  are 

n{t)  =Si{t-Tn)+S2it-'Cn)  (2.1) 

and 

riit)  =  Sl{t-T2l)+S2it-T:22),  (2.2) 

where  tmn  represents  the  delay  between  receiver  m  and  transmitter  n.  Any  differences  in 
signal  losses  across  receivers  are  assumed  to  be  negligible. 

Applying  the  Fourier  Transform  to  (2.1)  and  (2.2)  gives 

Riif)  =  5i(/)e-^“^'i  +S2(/)e-^“^'2  (2.3) 

and 

Riif)  =  5i(/)e-^'“"2i  (2.4) 
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where  co  =  2nf  and  S„(/)  and  Rm{f)  are  the  Fourier  transforms  of  the  transmitted  and 
reeeived  signals,  respeetively.  In  matrix  form,  and  extended  for  m  reeeivers  and  n  transmit¬ 
ters,  (2.3)  and  (2.4)  beeome 


Ti(/)  ■ 

.  RM)  , 

(2.5) 


Given  adequate  estimates  of  the  delays  Tm„,  (2.5)  ean  be  viewed  as  the  eanonieal  Ax  =  b 
linear  algebra  problem  for  all  points  in  /,  where  A  and  b  are  known,  and  linear  algebra 
teehniques  ean  be  applied  to  reeover  x,  whieh  holds  the  frequeney-domain  deseription  of 
eaeh  transmitted  signal. 


Working  with  sampled  data  of  length  N,  the  diserete  Fourier  transform  (DFT),  typieally  im¬ 
plemented  via  the  fast  Fourier  transform  (FFT),  yields  N  frequeney  bins,  with  the  frequeney 
of  eaeh  bin  given  by 

(O  =  2;r^/„  (2.6) 

in  radians  per  seeond  (rad/see),  where  i  is  the  bin  number  and  fs  is  the  sampling  frequeney, 
in  Hertz  (Hz).  This  means  that  there  are  N  separate  instanees  of  (2.5)  to  solve.  It  is  eon- 
venient  to  reshape  the  problem  from  N  unique  matrix  equations  to  a  single  equation  that 
ineorporates  all  N  frequeneies.  In  that  ease,  the  exponential  delay  terms  are  arranged  in 
mxn  bloeks  on  the  diagonal  of  a  sparse  matrix  and  (2.5)  becomes 


■  A(l) 

0 

■  s(i)  ■ 

■  R(l)  ■ 

0 

MN)  _ 

.  . 

(2.7) 


where 


7^^  jj  fs'^Xn 


A(i)=  ;  ;  (2.8) 

and  S  and  R  are  column  vectors  of  S(/)  =  [Si{i)  and  R(/)  =  [^i(/) 
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For  convenient  referenee,  (2.7)  is  summarized  as 


AS  =  R  (2.9) 

where  A,  written  without  an  index,  refers  to  the  eomplete  set  of  A(/)  blocks  arranged  on 
the  diagonal,  as  in  (2.7),  while  S  and  R  are  the  eomplete  eolumns  of  S  and  R,  respeetively. 


2.2  Geometry  and  Delay  Overview 


Lry{2)- 

Ltyi'2)  — 


S' 


3  3 


(a)  Top  view 


Figure  2.1.  Geometry  determines  the  signal  delay. 


The  delay  T  in  the  above  equations  is  the  propagation  delay,  whieh  is  a  funetion  of  the 
geometric  relationship  of  the  transmitters  and  reeeivers.  An  example  configuration  of  two 
transmitters  and  two  reeeivers  is  depieted  in  Figure  2.1,  but  this  is  readily  extended  to 
aeeommodate  additional  transmitters  and  reeeivers.  For  simplicity  of  the  distanee  ealeula- 
tion,  this  eonfiguration  assumes  that  the  receivers  travel  on  parallel  paths,  though  there  is 
no  requirement  for  sueh  a  geometry.  The  eonfiguration  measures  transmitter  and  reeeiver 
loeation  from  an  origin  point  and  reeords  those  locations  either  in  eolumn  vectors  Lrx  and 
Lry  for  receivers  or  Ltx  and  Lty  for  transmitters.  The  perpendieular  slant  distanee  between 
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the  receiver  path  and  transmitter  is  computed  first,  using  the  law  of  cosines  as 


dxmn  =  \l  {altm+Ref  +Rl-  2(a/t„  +7?e)7?eCOS  (0), 


(2.10) 


where 


e  = 


Lrxij^^  Ltxijd) 

~R. 


(2.11) 


is  the  interior  angle  between  the  receiver  path  and  the  transmitter  location,  as  depicted  in 
Figure  2.1b.  With  the  perpendicular  distance,  the  angle  7  is  defined  by 


tan  7  = 


Lf-yi/n)  Lfyin) 


(2.12) 


and  the  desired  distance  dmn  can  be  determined  using 


cos  7  = 


(2.13) 


Substituting  (2.12)  into  (2.13),  the  direct  distance  between  the  transmitter  and  receiver  is 
provided  by 


dmn  — 


dx 


COS 


tan 


—  \  (  (^)  (^^ 


(2.14) 


The  propagation  delay  is  then  calculated  using  Xrnn  =  dmnjc,  where  c  is  the  speed  of  light. 
A  MATLAB  implementation  of  this  approach  is  provided  in  Appendix  C.l. 


A  straightforward  solution  to  the  canonical  Ax  =  b  problem  requires  an  invertible  A  ma¬ 
trix.  In  this  scenario,  that  requires  both  a  cooperative  geometry  and  an  equal  number  of 
transmitters  and  receivers.  Geometry  can  be  evaluated  in  advance,  which  should  enable 
the  identification  and  avoidance  of  unfavorable  collection  geometries.  The  approach  is  first 
evaluated  in  the  m  =  n  case,  but  the  proliferation  of  transmitters  coupled  with  the  limited 
receiver  assets  makes  the  m  <  n  condition  worthy  of  additional  consideration. 


2.3  Narrowband  Evaluation 

For  m  =  n,  the  transmitted  signals  can  be  recovered  as  long  as  the  matrix  of  phase  delays 
is  invertible.  Such  a  geometry  was  used  to  evaluate  the  approach  using  narrowband  signals 
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in  an  effort  to  evaluate  the  validity  of  the  approaeh. 


2.3.1  Test  Geometry 

The  values  used  in  the  following  analysis  are  summarized  in  Table  2. 1 .  The  reeeiver  loea- 
tions  along  the  x-axis  were  ealeulated  to  give  a  look  angle  a  of  20°  from  nadir.  Locations 
along  the  y-axis  were  selected  to  avoid  complementary  delays.  While  such  delays  would 
have  resulted  in  a  delay  matrix  with  orthogonal  eigenvectors  and  simplified  the  analysis 
somewhat,  that  is  not  necessarily  representative  and  was,  therefore,  avoided. 

Table  2.1.  Summary  of  geometry  parameters. 


Description 

Parameter 

Value 

Number  of  transmitters 

n 

2 

Number  of  receivers 

m 

2 

Receiver  altitude 

alt 

400  km 

Earth  radius 

Re 

6371  km 

Look  angle  from  nadir 

a 

20° 

Transmitter  location  (x) 

Ltx 

[ 0  0 ] km 

Transmitter  location  (y) 

Rty 

[  0  30  ]  km 

Receiver  location  (x) 

Rrx 

[  146.2  146.2]  km 

Receiver  location  (y) 

Lry 

[  0  50]  km 

The  geometry  calculations,  as  described  in  Section  2.2,  and  construction  of  the  A  matrix 
defined  in  (2.8)  were  implemented  in  a  MATLAB  functions  delays_f  rom_laydown()  and 
block_deltas_f  rom_delays  () ,  which  are  provided  in  Appendix  C.l  and  Appendix  C.2. 

2.3.2  Signal  Generation  and  Recovery 

Once  the  delays  between  transmitter  and  receiver  pairs  Tmn  are  established,  it  is  possi¬ 
ble  to  generate  the  test  signals  as  both  transmitted  and  received.  The  MATLAB  function 
calc_signals_nb()  is  used  to  generate  both  the  transmitted  and  received  signals  and  is 
included  in  Appendix  C.3.  For  the  narrowband  tests,  signals  were  generated  of  the  form 

Snit)  =A„COs{(Ont)  (2.15) 

where  A„  is  the  signal  amplitude  for  the  nth  transmitter  (A„  =  1  for  all  of  the  cases  discussed 
below)  and  (0„  =  2k  fn  is  the  corresponding  signal  frequency.  For  the  first  test,  signals  were 
selected  that  consisted  solely  of  a  single  fundamental  frequency  of  the  Fourier  transform 
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domain  -  that  is,  frequencies  that  conform  to  integer  values  of  i  in  (2.6).  In  that  case, 
an  integer  number  of  waveforms  appears  in  the  time-domain  signal,  and  the  frequency- 
domain  representation  has  two  non-zero  values,  corresponding  to  That  restriction 

was  not  applied  to  tOi  in  the  second  narrowband  test  described  below.  Specific  parameters 
for  each  test  are  provided  in  Table  2.2.  The  transmitted  signals  were  then  sampled  at  Fs, 
and  a  corresponding  frequency-domain  set  of  samples  was  created  using  the  FFT. 

The  sampled  received  signals  can  be  constructed  as  the  sum  of  the  transmitted  signals, 
appropriately  delayed,  while  neglecting  any  attenuation,  and  ignoring  any  noise  effects 
temporarily.  As  the  transmitted  signals  are  defined  in  continuous  time,  the  delayed  versions 
can  be  computed  directly  for  each  of  the  sample  points.  Each  received  signal  is  then 

n 

'^mi)  1  (2.16) 

i=l 

and  the  corresponding  frequency-domain  signal  is  computed  with  the  FFT. 


Table  2.2.  Signal  and  sampling  parameters,  narrowband  tests. 


Description 

Parameter 

Value 

Common 

Sampling  frequency 

Fs 

1.024  MHz 

Number  of  samples 

N 

32 

Signal  magnitude 

1 

First  NB  Test 

Frequencies  (Hz) 

fl 

32,000 

fl 

187,500 

Second  NB  Test 

Frequencies  (Hz) 

fl 

47,685 

fl 

187,500 

Having  generated  received  signals  that  incorporate  the  properly  delayed  transmissions,  it  is 
possible  to  attempt  recovery  of  the  transmitted  signals  by  solving  (2.7)  for  S.  The  MATLAB 
left-division  operator  is  used  to  determine  that  solution  in  the  recover_f rom_delays  () 
function,  included  in  Appendix  C.4.  Solving  (2.7)  yields  the  frequency-domain  description 
of  the  transmitted  signals.  The  corresponding  time-domain  signals  are  then  computed  using 
the  inverse  fast  Fourier  transform  (IFFT). 
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2.3.3  Performance  Evaluation 

In  addition  to  a  qualitative  comparison  of  the  original  and  recovered  signals  in  both  the  time 
and  frequency  domains,  some  quantitative  assessment  of  the  recovered  signals  is  useful. 
Those  methods  are  described  here. 


Power  Angle 

As  the  test  signals  used  are  real- valued,  the  recovered  signals  should  share  that  property.  By 
comparing  the  signal  power  contained  in  the  real  versus  the  imaginary  part  of  the  recovered 
signal,  we  can  compute  a  power  angle,  the  orientation  of  the  signal  power  in  the  complex 
plane. 


For  a  sampled  signal,  the  power  is  given  by 


(2.17) 


where  N  is  the  number  of  samples  in  the  signal  s,  which  is  a  column  vector,  and  s*  is 
the  conjugate  transpose  of  the  column,  such  that  s*s  computes  the  inner  product  of  the 
signal  [19].  Signal  power  is  composed  of  real  and  imaginary  power.  The  angle  of  the 
power  vector  in  the  complex  plane  is  given  by 


&p  =  tan 


/3{s}^3{s}\ 


(2.18) 


where  the  real  and  imaginary  parts  of  the  signal  are  denoted  with  the  9l{}  and  3{}  opera¬ 
tors,  respectively,  and  the  length  factors  (l/N)  from  (2.17)  cancel.  A  &p  or  power  angle  (as 
referenced  below)  of  zero  indicates  that  the  signal  power  is  completely  real.  Non- zero  val¬ 
ues  indicate  that  the  recovery  process  resulted  in  a  signal  that  has  some  amount  of  reactive 
power,  which  is  unexpected,  and  serves  as  a  first  check  of  the  recovered  signal. 


Error  Figure 

An  understanding  of  the  differences  between  the  original  and  recovered  signal  is  also  a 
useful  indicator  of  performance.  Such  an  understanding  is  given  by  the  energy  or  power  of 
the  difference  between  the  two  signals,  particularly  when  normalized  by  the  energy  of  the 
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original  signal.  This  normalized  error  energy  or  error  figure  is  given  by 

Error  Figures  (2.19) 

s*s 

where  s  represents  the  reeovered  signal  and  s  the  original  signal  [19].  Zero  error  figure 
indieates  that  there  is  no  energy  in  (s  —  s)  —  that  is,  there  is  no  energy  in  the  reeovered 
signal  that  was  not  in  the  original  signal,  whieh  is  the  ideal  ease. 

Scalar  Product  or  Correlation  Coefficient 

The  similarity  of  two  signals  ean  be  assessed  by  eomputing  the  sealar  produet,  whieh  is 
maximal  when  the  signals  are  proportional  [12].  Sehwarz’s  inequality  relates  the  sealar 
produet  to  the  signal  powers,  sueh  that,  in  the  ease  of  transmitted  and  reeovered  signals,  s 
and  s, 

s*s  <  (2.20) 

where  the  power  of  the  diserete  time  signal  is  given  by  (2.17).  To  faeilitate  eomparison  of 
signals  with  different  power  levels,  it  is  desirable  to  normalize  to  unit  power  signals,  sueh 
that 

;t(s,s)  =  ^=.  (2.21) 

The  eorrelation  eoeffieient  ean  vary  from  -1  to  1,  where  a  value  of  1  indieates  proportional 
signals,  0  indieates  the  signals  are  orthogonal,  and  -1  indieates  that  the  signals  are  propor¬ 
tional  but  with  differing  sign  eonvention  [20]. 

Carlson  [12]  identifies  that  the  eross-eorrelation  funetion  is  the  sealar  produet  of  two  sig¬ 
nals  with  the  seeond  signal  delayed  by  t,  and  the  eross-eorrelation  funetion  measures  the 
similarity  of  the  two  signals  as  a  funetion  of  that  delay.  For  purposes  of  this  effort,  we 
are  interested  in  the  eorrelation  eoeffieient  when  T  =  0,  whieh  eorresponds  to  the  sealar 
produet.  Signals  are  deseribed  as  well  eorrelated  when  the  normalized  sealar  produet  ap- 
proaehes  one,  and  eorrelation  is  assumed  to  be  at  t  =  0. 

Subspace  Angle 

Taking  a  page  from  linear  algebra,  we  imagine  the  sampled  signals  as  veetors  in  N- 
dimensional  spaee.  Eaeh  veetor  deseribes  a  subspaee  eontained  within  that  A-dimensional 
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space.  In  that  context,  identical  signals  are  parallel  to  each  other,  pointing  in  the  same 
direction  in  that  A^-dimensional  space.  The  angle  between  two  such  vectors  is  given  by 


0  =  cos  ^ 


(2.22) 


where  ||  ■  ||  denotes  the  function  that  computes  the  length  or  2-norm  of  a  vector.  A  0  =  0 
result  indicates  parallel  vectors  while  a  0  =  90°  indicates  completely  orthogonal  vectors. 

Recognizing  that  Ps=  ||s|p/A  =  s*s/A  shows  that  subspace  angle  is  the  inverse  cosine  of 
the  normalized  scalar  product.  Throughout  this  document,  the  subspace  angle  is  reported 
in  degrees  as  an  alternative  perspective  to  the  scalar  product. 


2.3.4  Results 


Frequency (kHz) 


Frequency  (kHz) 


(a)  First  signal  (b)  Second  signal 

Figure  2.2.  Frequency-domain  plot  of  narrowband  signals  as  transmitted  and  recovered,  /i  =32 
kHz. 


Table  2.3.  Performance  figures  for  /i  =  32  kHz. 

Signal  1  Signal  2 

Power  Angle  (°)  5. 1252  x  5.3304  x 

Error  Figure  1.5111x10^26  i.i  12  x  10^26  ^  17^06  x 

Subspace  Angle  (°)  0  1.8149  x  lO^^^  _  1.8149  x 

Correlation  1  13-1.0033x  10^21/ 

The  results  of  the  first  narrowband  test  indicate  a  very  strong  agreement  between  the  trans¬ 
mitted  signals  and  the  numerically  recovered  transmitted  signals.  The  frequency-domain 
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plots^  (including  both  the  magnitude  and  phase  of  each)  of  the  transmitted  and  recovered 
signals  for  the  first  test,  with  f\  =  32  kHz,  are  provided  in  Figure  2.2.  For  clarity,  the 
phase  plots  only  include  data  points  for  those  frequencies  whose  magnitude  is  at  least  1% 
of  the  maximum  value  in  the  magnitude  plot.  Both  recovered  signals  demonstrate  good 
agreement  with  the  transmitted  signal,  including  very  small  phase  error  for  the  non-zero 
frequency  components  of  the  signal.  The  corresponding  performance  figures  for  this  test 
are  provided  in  Table  2.3.  The  power  angle  results  show  that  each  of  the  recovered  sig¬ 
nals  is  almost  entirely  real,  with  imaginary  parts  on  the  order  of  computer  rounding  error, 
Cmachine-  Similarly,  this  test  yields  very  small  error  figures,  indicating  that  there  is  very  lim¬ 
ited  power  difference  between  the  two  signals.  Lastly,  the  angle  between  these  two  vectors, 
when  imagined  in  A^-dimensional  space,  is  very  nearly  zero,  indicating  parallel  vectors, 
which  is  the  ideal  case.  Overall,  this  indicates  that  the  technique  has  some  promise,  and  a 
good  correlation  coefficient  in  the  time  domain  is  expected. 


(a)  First  test  (b)  Second  test 

Figure  2.3.  Narrowband  tests,  transmitted  and  recovered  signals. 

The  time-domain  results  for  the  first  narrowband  test  are  provided  in  Figure  2.3a.  As 
expected,  the  recovered  signal  is  closely  matched  to  the  original  signal.  Similarly,  the 
time-domain  signal  results  from  the  second  test,  using  f\  =  47,685  Hz  are  provided  in 
Figure  2.3b.  That  case  shows  considerable  disagreement  between  the  original  and  recov¬ 
ered  signals.  Examination  of  the  frequency-domain  plots,  provided  in  Figure  2.4,  provide 
insight  into  the  problem. 

^All  signals  throughout  the  document  are  simulated  and  assumed  to  be  either  voltage  or  current.  Units  are 
left  unspecified  in  all  figures. 
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Frequency (kHz)  Frequency (kHz) 


(a)  First  signal  (b)  Second  signal 

Figure  2.4.  Frequency-domain  plot  of  narrowband  signals  as  transmitted  and  recovered,  /i  = 
47,685  Hz. 


Table  2.4.  Performance  figures  for  /i  =47,685  Hz. 


Signal  1 

Power  Angle  (°)  0.00213 

Error  Figure  10.34 

Subspaee  Angle  (°)  110.5  —  0.01644/ 

Correlation  -0.3497  4-0.0002688/ 


Signal  2 

0.0003375 

6.676 

72.01 -F 2.3432  xl0-l^/ 
0.3088-3.8898xl0-l‘^/ 


Unlike  the  first  test,  in  the  seeond  test,  f\  is  transmitted  as  a  sinusoid  that  maps  onto  many 
DFT  indiees,  sinee  f\  ^  n/Tsamp  (where  is  n  is  an  integer),  and  this  transmitted  signal 
demonstrates  a  smooth  phase  throughout  most  of  the  speetrum.  (Strum  refers  to  this  as 
leakage  [21].)  The  reeovered  signal  laeks  this  smooth  phase,  instead  yielding  a  mueh  more 
erratie  phase.  While  the  primary  frequeney  of  the  seeond  signal  appears  in  the  reeovered 
version  (though  with  redueed  magnitude),  there  are  a  number  of  other  frequeney  eompo- 
nents  present  in  the  reeovered  version  that  are  not  present  in  the  transmitted  signal.  The 
performanee  figures,  provided  in  Table  2.4,  also  indieate  poor  reeovery  performanee.  While 
the  power  angle  remains  low  (indieating  real  valued  signals),  the  error  figure,  subspaee  an¬ 
gle,  and  eorrelation  all  indieate  little  agreement  between  the  original  and  reeovered  signals. 


2.4  Frequency  Response  Analysis 

The  approaeh  outlined  in  Seetion  2.1  depends  upon  the  faet  that  a  known  delay  in  the  time 
domain  translates  to  a  predietable  phase  shift  in  the  frequeney  domain  and  that  the  Fourier 
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Transform  allows  for  transition  between  those  domains.  While  these  are  known  properties 
of  the  Fourier  Transform,  the  results  described  above  indicate  that  there  are  likely  some 
details  of  those  properties  that  were  not  adequately  considered. 


2.4.1  Frequency  Response  Definition 


In  order  to  study  the  necessary  relationship,  the  delay  was  considered  as  a  frequency  re¬ 
sponse,  such  that 


H{0)) 


(2.23) 


where  ,^[-]  represents  the  Fourier  transform  operation.  By  computing  the  FFT  of  a  signal 
and  a  delayed  copy,  it  is  possible  to  compare  the  resulting  frequency  response  with  the 
expected  result.  As  the  approach  outlined  in  Section  2.1  is  focused  on  the  phase 

relationships,  the  frequency  response  magnitude  is  less  important.  Any  deviations  from 
unit  magnitude  are  canceled  by  the  IFFT  operation.  The  phase  relationship  is  the  critical 
component,  as  without  a  linear  phase  shift  over  varying  frequency,  the  approach  described 
in  Section  2.1  breaks  down  quickly.  For  all  of  the  tests  below,  the  delay  is  T  =  IT^amp, 
where  Tsamp  is  the  sampling  interval. 


2.4.2  Fourier  Transform  Relationships 

The  delay  property  applied  in  (2.3)  and  (2.4)  is  a  property  of  the  continuous  Fourier  trans¬ 
form  (CFT).  For  sampled  data,  a  similar  property  exists  for  the  discrete  time  Fourier  trans¬ 
form  (DTFT).  When  the  sampled  data  is  A-periodic,  the  DTFT  can  be  simplified  to  the 
DFT,  which  the  FFT  implements  efficiently  and  also  shares  a  shifting  property  that  yields  a 
linear  phase  shift  with  frequency.  The  important  characteristic  of  the  DFT,  however,  is  that 
the  original  samples  must  be  A-periodic  in  order  for  it  to  be  accurately  applied  in  place  of 
the  DTFT.  Ultimately,  this  requirement  for  A-periodic  data  means  that  the  shifting  property 
of  the  DFT  is  a  circular  shift  [21].  For  data  that  is  not  A-periodic,  this  difference  yields  a 
different  phase  response  for  the  frequency  response,  and  without  somehow  accounting  for 
this  difference,  the  above  approach  breaks  down.  These  differences  are  illustrated  below. 


2.4.3  Narrowband  Cases 

A  narrowband  signal  of  the  form  described  by  (2. 15)  was  generated  with  32  sampled  points. 
The  magnitude  and  phase  of  the  frequency  response  are  provided  in  Figure  2.5  along  with 
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a  depiction  of  the  signals  with  two  different  delay  types.  The  first  delay,  labeled  as  “Con¬ 
tinuous  time  shift”  was  generated  by  sampling  a  delayed  continuous  time  signal,  while  the 
second  delay,  labeled  as  “Circular  shift”  was  computed  with  a  circular  shift  of  the  original 
signal.  Both  signals  were  delayed  by  two  sample  intervals.  Frequency  responses  for  both 
delay  types  as  well  as  the  ideal  response  are  provided  in  Figure  2.5a.  Note  that 

while  not  required  in  this  case,  later  cases  involve  magnitude  deviations  from  one  on  the 
order  of  machine  error.  For  consistency,  the  magnitude  plot  is  plotted  as  |//((0)|  —  1  to 
clearly  depict  those  small  deviations. 


(a)  Frequency  response  (b)  Delay  comparison 

Figure  2.5.  Delay  analysis  of  the  /i  =  32  kHz  signal. 


Recall  from  Section  2.3.4  that  this  signal  provided  excellent  recovery  performance.  While 
the  phase  plot  does  not  show  good  correlation  between  either  of  the  frequency  responses 
and  the  ideal  response,  it  is  important  to  recognize  that  in  the  frequency  domain,  ideally 
this  signal  has  zero  magnitude  for  all  but  two  frequencies.  When  the  frequency  response 
is  computed,  both  the  dividend  and  divisor  are  approaching  £machine  for  all  but  these  two 
frequencies,  and  the  result  shows  considerable  variation.  The  two  frequencies  with  non¬ 
zero  magnitude  are  marked  with  small  squares,  and  each  of  these  data  points  show  good 
agreement  across  the  two  delay  types  and  the  expected  frequency  response.  The  time- 
domain  view  of  the  two  delay  types  shows  they  yield  similar  results,  as  at  this  frequency 
and  sample  rate,  the  signal  is  //-periodic. 

The  second  narrowband  case,  as  shown  previously,  has  a  markedly  poorer  recovery  perfor- 
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(a)  Frequency  response  (b)  Delay  comparison 

Figure  2.6.  Delay  analysis  of  the  /i  =47.7  kHz  signal. 


mance,  and  the  frequency  response  analysis  of  that  case,  provided  in  Figure  2.6,  illustrates 
the  problem  underlying  that  performance.  As  the  time-domain  signal  in  this  case  is  not 
composed  of  a  single  component  in  the  discrete  frequency  domain,  the  magnitude  of  the 
various  frequencies  is  non-zero  and,  therefore,  yields  a  much  less  erratic  frequency  re¬ 
sponse  for  both  types  of  shifting.  While  the  circular  shift  of  the  signal  does  appear  well 
correlated  with  the  ideal  case,  the  signal  that  was  shifted  in  the  continuous  time  domain 
prior  to  sampling  has  a  radically  different  frequency  response,  particularly  with  respect 
to  the  phase  angle.  The  difference  results  from  the  fact  that  the  signal  is  not  A-periodic, 
which  is  clear  from  the  difference  in  the  shifted  signals  presented  in  Figure  2.6b.  There, 
two  samples  (corresponding  to  the  T  =  ITsamp  delay)  of  the  discrete  (or  circular)  shift  are 
not  correlated  with  the  corresponding  samples  of  the  continuous  time  shifted  signal.  The 
delayed  signals  (as  received)  were  generated  using  a  time  shift  of  a  function  defined  in 
continuous  time.  A  correct  solution  for  S  from  (2.7)  depends  upon  a  phase  shift  that  would 
have  resulted  from  a  circular  shift  vice  a  delay  of  the  continuous  time  signal.  An  accurate 
recovery  requires  that  the  effects  of  the  continuous  time  delay  closely  resemble  those  of  a 
circular  shift. 

2.4.4  Wideband  Cases 

For  comparison  to  the  narrowband  case,  a  wideband  test  case  was  also  conducted,  using 
an  additive  white  Gaussian  noise  (AWGN)  signal.  To  facilitate  a  more  direct  comparison 
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between  this  and  the  narrowband  case,  this  signal  was  also  generated  using  32  sampled 
values.  As  AWGN  is  not  a  time-dependent  signal,  the  continuous  time  delay  was  created  by 
prepending  two  samples  of  AWGN  to  the  front  of  the  signal,  while  the  discrete  time  delay 
was  computed  using  a  circular  shift,  as  in  the  narrowband  case  above.  Using  the  same 
format  as  above,  we  provide  transfer  functions  for  both  delay  types  as  well  as  the  ideal 
response  in  Figure  2.7a.  The  time-domain  plot  of  the  three  signals  is  provided  in 

Figure  2.7b,  and  as  in  the  second  narrowband  case,  also  exhibits  two  points  in  the  circularly 
shifted  signal  that  do  not  match  those  in  the  continuous  time  delayed  signal.  Despite  that 
difference,  the  phase  plot  of  Figure  2.7a  shows  that  while  the  continuous  time  shift  line  does 
exhibit  a  few  points  of  departure  from  the  ideal  and  circular  shift  cases,  there  is  generally 
much  better  agreement  between  the  continuous  time  shift  curve  and  the  ideal  case  than  in 
the  narrowband  cases,  with  the  overall  trend  of  the  curve  following  that  of  the  ideal  curve. 
This  indicates  that  a  wider  bandwidth  signal  might  adequately  support  the  signal  recovery 
approach  described  above. 


(a)  Frequency  response  (b)  Delay  comparison 

Figure  2.7.  Delay  analysis  of  the  AWGN  signal,  N  =  32. 


The  same  test  was  repeated  with  additional  samples  to  determine  if  the  performance  was 
a  function  of  the  number  of  samples  of  the  AWGN  signal.  With  additional  samples,  the 
number  of  samples  that  differ  between  the  continuous  delay  and  circular  shift  is  a  smaller 
fraction  of  the  overall  signal.  An  AWGN  signal  was  generated  using  N  =  2048,  and  the 
plots  of  the  frequency  response  performance  are  provided  in  Figure  2.8a  while  Figure  2.8b 
depicts  the  first  32  samples  of  the  complete  time-domain  signal  in  order  to  more  clearly 
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show  the  two  different  points  and  the  overall  two-sample  delay  of  the  signals.  While  there  is 
some  phase  noise  present  in  the  eontinuous  time-shifted  signal,  it  is  important  to  recognize 
that  this  signal  is  plotted  over  the  [— ;r, -l-;r]  interval  which,  in  some  cases,  causes  what 
might  otherwise  be  a  smaller  error  to  appear  more  significant.  Aside  from  this  artifact,  the 
continuous  time  plot  does  appear  more  closely  aligned  with  the  ideal  case  than  any  of  the 
previous  frequency  response  tests.  This  improvement  is  likely  due  to  the  larger  number  of 
samples.  With  more  samples,  the  two  points  that  are  impacted  by  the  shift  differences  are  a 
much  smaller  fraction  of  the  total  signal  and  have  little  impact  on  the  frequency  response. 


(a)  Frequency  response  (b)  Delay  comparison 

Figure  2.8.  Delay  analysis  of  the  AWGN  signal,  N  =  2048. 


2.4.5  Zero  Padding 

To  further  improve  the  similarity  between  the  delay  that  occurs  in  continuous  time  and  the 
circular  shift  that  is  required  to  realize  a  phase  shift  in  the  frequency  domain  of  the 
sampled  signal,  zero  padding  and  windowing  can  be  applied  to  the  collected  signal. 

Intuitively,  if  a  signal  is  transmitted  for  a  limited  period  of  time,  a  receiver  receives  no  signal 
while  the  leading  edge  of  the  signal  is  propagating,  and  the  signal  arrives  only  after  the 
propagation  delay  T.  The  causal  nature  of  the  linear  time  invariant  (LTI)  channel  requires 
that  the  samples  of  the  signal  arrive  in  the  same  order  in  which  they  were  transmitted. 
Without  zero  padding,  the  circular  shift  associated  with  that  delay  in  the  discrete  time 
domain  results  in  some  fraction  of  the  samples  from  the  trailing  edge  of  the  signal  shifting 
to  the  leading  edge,  which  is  not  representative  of  the  actual  signal  received  after  the  delay. 
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By  padding  a  time-limited  signal  with  trailing  zeros,  the  propagation  delay  results  in  some 
fraetion  of  those  zeros  moving  to  the  start  of  the  signal  rather  than  some  portion  of  the 
trailing  edge  of  the  original  signal,  preserving  the  eausal  nature  of  the  transmission  and 
propagation  delay. 

Mathematieally,  for  a  sampled  signal 

x[n]  =x{nT)  forn  G  [0,A/'—  1]  (2.24) 


of  duration  N  and  integer  sample  numbers  n,  zero  padding  yields  a  signal 


x[n\  =  x{nT)I^Qi^^l-^{n)  forne  [0,A^— 1], 


(2.25) 


where 


1 ,  if  n  G  A 
0,  otherwise 


(2.26) 


is  known  as  the  indieator  funetion  and  N  is  the  new  length  of  the  signal.  Adequate  zero 
padding  requires 

(2.27) 

where  T  is  the  duration  of  the  signal,  Tj  is  the  propagation  delay  and  Tsamp  is  the  sampling 
interval. 


The  proeeeding  evaluation  of  shifting  types  was  repeated  with  zero  padded  AWGN  signals 
for  signal  lengths  of  A  =  32  and  N  =  2048.  The  signal  delay  evaluated  in  this  ease  is 
ideal,  sueh  that  it  is  an  integer  number  of  samples,  the  transmitted  signal  is  time-limited, 
and  the  delayed  signal  ineorporates  the  eomplete  transmitted  signal.  The  eomparison  of 
the  eontinuous  time  shift  and  eireular  shift  are  provided  in  Figure  2.9  for  the  N  =  32  ease 
and  Figure  2.10  for  the  N  =  2048  ease.  Regardless  of  signal  length,  the  eireular  shift  and 
eontinuous  time  delay  are  nearly  identieal,  whieh  demonstrates  that  under  ideal  eonditions, 
the  eontinuous  time  delay  that  results  from  signal  propagation  ean  be  removed  through 
either  a  eireular  shift  in  the  diserete  time  domain  or  by  multiplieation  of  a  phase  shift  in  the 
frequeney  domain. 
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(a)  Frequency  response  (b)  Delay  comparison 

Figure  2.9.  Delay  analysis  of  the  AWGN  signal,  N  =  32,  with  zero-padding. 


(a)  Frequency  response  (b)  Delay  comparison 

Figure  2.10.  Delay  analysis  of  the  AWGN  signal,  N  =  2048,  with  zero-padding. 


2.4.6  Phase  Error  Calculation 


Rather  than  rely  simply  upon  a  qualitative  assessment,  it  is  preferable  to  evaluate  the  phase 
errors  deseribed  above  quantitatively.  For  both  the  continuously  shifted  and  circularly 
shifted  signals,  the  normalized  phase  error  is  given  by 


Phase  Error  = 


ZH{co)-ZH{co)\\ 


(2.28) 
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where  ZH{co)  is  the  phase  angle  of  the  frequeney  response  for  the  shifted  signal  (either 
eontinuous  time  or  eireularly  shifted)  and  ZH{(o)  is  the  phase  angle  of  the  ideal  transfer 
funetion,  =  —cot.  Ideally,  this  phase  error  is  zero.  The  phase  angle  results  are 

presented  in  Table  2.5.  Reeall  that  in  the  narrowband  ease  with  /  =  32  kHz,  both  the 
delay  approaehes  had  eonsiderable  error  as  a  result  of  the  numerous  frequeneies  with  near 
zero  response.  Other  than  that  ease,  the  error  for  the  eireular  shift  approaeh  is  roughly 
Bmachine^  wMoh  eonfirms  the  graphieal  assessment  above.  For  the  eontinuous  time  delayed 
signals,  in  the  narrowband  oases,  the  phase  error  is  approximately  one,  whieh  indioates 
that  the  phase  error  is  of  the  same  magnitude  as  the  ideal  phase  response  itself.  Some 
improvement  is  provided  in  the  wideband,  N  =  32  ease,  while  the  wideband,  N  =  2048  ease 
shows  a  eonsiderable  improvement,  validating  the  use  of  additional  points  with  an  AWGN 
signal  to  mitigate  the  differenees  between  the  eontinuous  and  diserete  time  shifts.  Onee  the 
signals  are  zero  padded,  the  eontinuous  time  delay  preeisely  mirrors  the  performanee  of 
the  eireular  shift.  It  is  important  to  reeognize  that  the  signals  used  in  the  zero  padded  eases 
were  time  limited  at  the  transmitter,  sueh  that  the  initial  samples  of  the  reeeived  signals  are 
zeros  eorresponding  to  the  propagation  delay  and  not  samples  of  the  signal  that  had  been 
transmitted  before  t  =  0  at  the  transmitter. 


Table  2.5.  Frequency  response  test  normalized  phase  errors. 

eireular  Continuous 


Narrowband,  /  =  32  kHz 
Narrowband,  /  =  47.7  kHz 
Wideband,  N  =  32 
Wideband,  N  =  2048 
Wideband,  N  =  32,  zero-padded 
Wideband,  N  =  2048,  zero-padded 


0.04727 

1.054 

1.5254  X  10-16 

1.017 

3.4068  X  10-1^ 

0.1097 

7.6214  X  10-1^ 

0.007304 

1.7453  X  10-16 

1.7453  X  10 

1.7632  X  10-16 

1.7632  X  10 

2.5  Wideband  Test 

Having  identified  a  signal  that  possesses  the  neeessary  phase  relationship  between  the  orig¬ 
inal  signal  and  one  that  has  been  shifted  in  eontinuous  time,  we  ean  evaluate  the  ability  to 
reeover  that  signal.  For  eontinuity,  the  test  geometry  remains  as  deseribed  in  Seetion  2.3.1. 
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2.5.1  Signal  Generation  and  Recovery 


Sampled  baseband  AWGN  signals  were  generated  for  each  transmitter.  The  signals  arrive 
at  the  different  receivers  with  delays  dependent  upon  the  transmitter  and  geometry.  As 
those  delays  are  based  on  the  specific  geometry  of  the  transmitters  and  receivers,  they  are 
not  an  integer  multiple  of  the  sampling  interval,  and  a  straightforward  circular  shift  of  the 
discrete  time  signals  is  inadequate  to  incorporate  this  precise  delay.  Two  methods  were 
considered  to  incorporate  the  appropriate  delay  into  the  transmitted  signals  —  delay  via 
phase  shift  in  the  frequency  domain  and  interpolation  of  the  sampled  time-domain  AWGN 
signals.  A  comparison  of  those  methods  is  provided  in  Appendix  A,  and  based  upon  that 
evaluation,  the  delay  was  applied  via  phase  shift  in  the  frequency  domain  of  the  transmitted 
signal  given  the  small  error  and  low  computational  complexity. 

As  discussed  in  Section  2.4.5,  it  was  necessary  to  incorporate  zero  padding  in  the  signals 
prior  to  computation  of  the  FFT.  The  minimum  padded  signal  length  is  given  by  (2.27), 
with  Tj  =  max(Tmn).  In  the  interest  of  FFT  computational  speed,  the  signal  was  padded  to 
the  next  largest  power  of  two,  such  that 

=  (2.29) 

where  [ ■]  denotes  the  ceiling,  or  round  up  to  the  nearest  integer,  operation. 

Given  the  unique  copies  of  the  signals  with  the  appropriate  delays  for  each  transmitter/re¬ 
ceiver  pair,  the  complete  received  signal  is  constructed  as 

rj^n]  =  £  |5[fc]e-A(2;r/yv)(Ar™) |  (2.30) 

i=\ 

where  ^  is  the  number  of  transmitters  and  n  is  the  time-domain  sample  index  number. 

The  MATLAB  function  calc_signals_wb()  used  for  signal  generation  is  provided  in 
Appendix  C.5,  while  the  basic  signal  parameters  used  are  provided  in  Table  2.6.  It  is 
important  to  note  that  the  received  signals  constructed  by  calc_signals_wb  are  delayed 
copies  of  a  time-limited  transmission,  as  opposed  to  signals  of  infinite  duration  that  are 


26 


sampled  for  a  specific  period  in  the  receiver.^ 

Table  2.6.  Signal  and  sampling  parameters,  wideband  tests. 

Description  Parameter  Value 

Sampling  frequency  Fs  1  MHz 

Number  of  signal  samples  N  2048 

Padded  signal  length  N  4096 

Desired  power  spectral  density  1  W  /  Hz 

Recovery  of  the  transmitted  signals  from  those  received  is  again  accomplished  by  solving 
(2.7)  for  S. 

2.5.2  Results 

Given  the  large  number  of  points  in  the  transmitted  signals,  plots  similar  to  those  used  in  the 
narrowband  test  are  less  informative  when  assessing  signal  recovery  performance.  Instead, 
subsets  of  the  time  and  frequency  domains  are  provided  in  Figure  2.11  and  Figure  2.12. 
These  show  excellent  correlation  between  the  two  signals,  with  the  exception  of  the  /  =  0 
points  in  both  Figure  2.12a  and  Figure  2.12b.  Those  points  are  not  recovered  when  solving 
(2.7)  as  the  equation  is  singular  at  A(l),  which  corresponds  to  0  Hz.^  (See  (2.8).)  Error 
plots,  depicting  the  difference  between  the  original  and  recovered  signals  are  presented  in 
Figure  2.13.  The  recovery  error  in  the  frequency  domain  is 

Ss=S-S,  (2.31) 

where  S  is  the  vector  of  frequency-domain  samples  of  the  original  signal  and  S  is  a  vector  of 
the  frequency-domain  samples  obtained  by  the  recovery  process.  Note  that  in  the  frequency 
domain,  the  error  at  /  =  0  has  been  omitted  from  the  plot  due  to  the  known  discrepancy  of 
the  recovery  process.  The  recovery  error  in  the  time  domain  is 

S5  =  s  — s,  (2.32) 

where  s  and  s  are  the  vectors  of  the  time-domain  samples  of  the  original  and  recovered 
signals,  respectively. 

^See  Section  5.2  for  further  discussion. 

^In  the  case  of  both  WiMax  and  LTE,  this  point  corresponds  to  the  DC  subcarrier,  which  both  standards 
specify  should  be  null.  The  recovery  process  simply  sets  this  point  to  zero  [7],  [22]. 
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The  performance  metrics  identified  in  Section  2.3.3  are  provided  in  Table  2.7.  Note  the 
small  Power  Angle,  Error  Figure  and  Subspace  Angle  figures,  and  the  correlation  coeffi¬ 
cient  very  near  one.  These  all  indicate  good  recovery  performance. 


Table  2.7.  Recovery  performance  figures  for  AWGN,  N  =  2048. 


Power  Angle  (°) 
Error  Figure 
Subspace  Angle  (°) 
Correlation 


Signal  1  Signal  2 

2.0427  X  10-^^  2.073  x  10~^^ 

6.0921  X  10^5  0.0001202-1.0917  x  10^44/ 

0.4472-4.263  X  10-12/  0.6283-1.6754x  10-12/ 

1  +  5.8073  X  10-16/  0.9999 -f  3.2064  X  10-16/ 
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Figure  2.11.  Wideband  test  signal,  time  domain. 


2.6  Performance  in  the  Presence  of  Noise 

To  be  a  useful  technique,  the  process  described  in  Section  2.1  must  demonstrate  reasonable 
performance  in  the  presence  of  noise.  To  evaluate  the  performance  with  noise,  the  above 
wideband  signal  test  using  AWGN  signals  was  repeated  with  varying  levels  of  AWGN 
noise  at  the  receivers.  No  filtering  was  applied  to  either  the  AWGN  signal  or  the  AWGN 
noise. 

As  above,  received  signals  were  constructed  using  (2.30)  with  determined  by  the  same 
transmitter  and  receiver  geometry  described  in  Section  2.3.1.  Noise  signals  were  con¬ 
structed  for  each  receiver  as  real  valued  random  samples  of  a  Gaussian  distribution  and 
scaled  by  the  standard  deviation  of  the  samples  in  order  to  ensure  unit  variance  of  the 
sampled  noise  signals. 
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Frequency  (kHz) 


Frequency  (kHz) 


Frequency  (kHz)  Frequency  (kHz) 


(a)  First  signal  (b)  Second  signal 

Figure  2.12.  Frequency  domain  of  wideband  test  transmitted  and  recovered  signals. 


(a)  Frequency  domain  (b)  Time  domain 

Figure  2.13.  Error  of  recovered  signals. 


At  each  signal-to-noise  ratio  (SNR)  test  point,  the  necessary  noise  power  level  was  com¬ 
puted  as 


Pr 

IQSNR/10’ 


(2.33) 


where  Pr  is  the  power  in  the  received  signal  and  the  desired  SNR  level  is  in  decibels  (dB). 
The  received  signal,  at  the  desired  SNR,  is 


^ m  —  “F 


(2.34) 
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The  signal  recovery  process  described  above  was  performed  and  the  correlation  coefficient 
between  the  original  and  recovered  signals  ;t^(s,s)  plotted  vs  SNR  in  Figure  2.14.  Some 
performance  degradation  is  expected,  as  the  recovery  process  solves  AS  =  R  for  the  signals 
S,  while  with  noise,  A5  +  T]  =  R  is  a  more  accurate  model  of  the  received  signals.  The 
approach  described  in  Section  2.1  makes  no  provision  for  identifying  or  eliminating  noise, 
and,  therefore,  even  under  ideal  performance  conditions,  the  recovered  signals  are  noisy 
versions  of  the  original  signals.  To  facilitate  an  understanding  of  how  the  recovery  process 
compares  to  what  might  be  considered  the  ideal  case.  Figure  2.14  includes  curves  for  the 
correlation  coefficient  of  the  original  signal  with  a  noisy  version  of  that  signal  ;t^(s,s), 
annotated  in  the  legend  as  “Original”,  where 

~  T  O'fira-  (2.35) 

As  this  is  intended  only  to  show  the  impact  of  noise  on  the  correlation  coefficient,  there 
was  no  attempt  to  include  this  noise  in  the  propagation  model. 

This  data  indicates  that  for  comparable  signal  correlation  coefficients,  the  recovery  process 
requires  roughly  35  dB  higher  SNR  than  is  necessary  if  the  recovery  process  is  not  uti¬ 
lized.  This  is  a  considerable  performance  degradation,  and  without  mitigation,  may  limit 
the  utility  of  this  approach.  See  Section  5.2  for  a  discussion  of  means  to  both  assess  this 
performance  impact  and  potential  means  to  mitigate  this  degradation.  Prior  to  expending 
the  effort  to  improve  performance  in  the  presence  of  noise,  the  recovery  method  should  be 
tested  with  the  OFDM  signals  used  by  fourth  generation  communication  systems. 


Figure  2.14.  Recovery  performance  versus  noise  level. 
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CHAPTER  3: 

Recovery  Using  Orthogonal  Frequency-Division 
Multiplexing  Signals 


In  order  to  be  useful  as  a  eo-ehannel  interferenee  mitigation  approaeh,  it  is  neeessary  to 
ensure  that  the  proeess  outlined  above  provides  adequate  performanee  with  representative 
signals.  OFDM  is  eentral  to  both  LTE  and  WiMax  teehnologies  whieh  offer  promise  as 
illuminators  for  a  passive  SAR  imaging  proeess.  A  brief  introduetion  to  the  OFDM  signal 
and  an  assessment  of  reeovery  performanee  with  OFDM  signals  is  presented  below. 


3.1  OFDM  Signal  Generation 

Fundamentally,  OFDM  ean  be  eonsidered  an  extension  of  frequeney-division  multiplex¬ 
ing  (FDM).  While  FDM  suffers  from  the  potential  for  interferenee  between  earriers  on 
differenee  frequeneies  (due  to  non-ideal  filters),  OFDM  avoids  this  with  the  use  of  orthog¬ 
onal  frequeneies,  sueh  that  for  a  symbol  duration  T  and  subearrier  frequeneies  fk  [12] 

I  cos(27t fkt)  cos{27t fjt)dt  =  0  Vy,  A:  where  7  y  A:.  (3.1) 

Given  the  trigonometrie  identity 

cos  6  cos  ^  ^  [eos(0  —  ^)  -|-eos(0  -I-  ^)] ,  (3.2) 

(3.1)  ean  be  rewritten  as 

^  cos[2m{fk-fj)]+cos[2m{fk  +  fj)]dt  =  0  Wj,k  where  j  ^  k.  (3.3) 

By  seleeting  fk+i  —  fk  =  Af  =  l/T,the  fk  —  fj  and  fk  +  fj  terms  are  an  integer  multiple  of 
A/,  the  integrals  are  zero  over  the  [0,  T]  interval  and  the  frequeneies  are  orthogonal. 

The  signal  that  represents  a  single  OFDM  symbol  is  assembled  from  a  series  of  K  subear- 
riers.  Eaeh  subearrier  is  a  sequenee  of  symbols  to  be  transmitted,  either  from  independent 
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sources  or  a  selection  of  K  symbols  to  be  transmitted  simultaneously  from  a  single  stream 
of  data.  Each  subcarrier  is  independently  modulated  (e.g.,  by  using  quadrature  phase-shift 
keying  (QPSK))  to  produce  K  streams  of  complex  data  Xi^[n].  For  a  given  OFDM  symbol, 
a  data  element  from  each  subcarrier  is  further  modulated  at  an  orthogonal  frequency  fk  and 
combined  with  the  other  modulated  subcarrier  signals,  as  in 

K-l 

w[n]  =  £  n  =  0,l,...K-\.  (3.4) 

k=0 

The  inverse  discrete  Fourier  transform  (IDFT)  provides  a  convenient  and  elegant  means  to 
implement  this  modulation  scheme  (see  [12]  for  a  detailed  discussion),  such  that 

w[n]=lD¥T{Xk[n]}.  (3.5) 

The  total  OFDM  signal  bandwidth  is 


B  =  KAf. 


(3.6) 


3.2  Recovery  Performance  with  OFDM 

Baseband  timelimited  OFDM  signals  were  generated  using  code  adapted  from  the  work  by 
Holder  [18].  The  parameters  for  each  set  of  test  signals  are  provided  in  Table  3.1.  (Note 
that  the  bandwidth  listed  there  is  a  two-sided  bandwidth.)  Those  signals  were  delayed 
and  combined  as  discussed  in  Section  2.5.1,  using  the  transmitter  and  receiver  geometry 
detailed  in  Table  2.1.  For  consistency  with  the  preceding  discussions,  the  OFDM  signals 
were  interpolated  and  resampled  at  1.024  MHz.  Relatively  low  bandwidth  OFDM  signals 
(as  compared  to  FTE  and  WiMax)  were  used  for  testing  to  keep  the  number  of  samples 
required  for  the  signal  delay  reasonable  at  this  sampling  rate.  See  Section  5.2  for  alternative 
approaches. 

Recovery  of  the  transmitted  signals  was  performed  by  solving  (2.7)  for  S.  Results  of  the 
recovery  process  are  provided  in  Table  3.2  through  Table  3.5.  Note  that  the  OFDM  signals 
are  not  generally  real- valued,  and  a  non-zero  power  angle  is  not  indicative  of  a  recovery 
shortfall.  A  power  angle  near  45°  indicates  comparable  real  and  imaginary  parts  of  the 
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Table  3.1.  OFDM  test  signal  parameter  summary. 


Set  1 

Set  2 

Set  3 

Set  4 

Subcarriers 

6 

128 

128 

128 

Modulation  Scheme 

QPSK 

QPSK 

QPSK 

BPSK 

Bandwidth  (kHz) 

600 

600 

1000 

1000 

OEDM  Symbols 

16 

4 

4 

10 

signal,  which  is  expected  with  the  OFDM  signals.  Subsets  of  the  signals  from  the  first 
parameter  set  depicted  in  the  time  and  frequency  domains  are  provided  in  Figure  3.1  and 
Figure  3.2.  Qualitatively,  these  indicate  good  recovery  performance  with  the  recovered  sig¬ 
nal  plotted  at  the  same  points  as  the  original  signal.  The  singularity  at  /  =  0  Hz  is  apparent 
in  the  frequency-domain  plots.  Time-  and  frequency-domain  plots  of  recovery  error  from 
the  first  set  of  signals  are  provided  in  Figure  3.3.  The  frequency-domain  error  indicates 
recovery  errors  on  the  order  of  machine  error,  with  a  few  exceptional  frequencies  that  have 
slightly  more  error  than  most.  The  steady-state  error  in  the  time  domain  is  likely  the  result 
of  the  singularity  at  /  =  0  Hz.  Each  set  of  OFDM  parameters  demonstrated  comparable 
performance,  with  some  minor  differences  that  are  likely  the  result  of  the  random  signal 
content  rather  than  a  function  of  one  of  the  OFDM  parameters.  The  complete  set  of  OFDM 
recovery  performance  plots  are  provided  in  Section  B.l.  This  consistent  and  accurate  re¬ 
covery  performance  suggests  that  the  method  is  likely  functional  across  a  large  collection 
of  OFDM-based  signals,  including  LTE  and  WiMax  transmissions.  It  is  important  to  rec¬ 
ognize  that  this  performance  was  examined  in  a  single  arrangement  of  transmitters  and 
receivers.  As  it  is  useful  to  understand  the  potential  impact  of  collection  geometry  on 
recovery  performance,  this  is  explored  further  in  the  following  chapter. 


Table  3.2.  Recovery  performance  figures  for  OFDM,  set  1. 


Power  Angle  (°) 
Error  Eigure 
Subspace  Angle  (°) 
Correlation 


Signal  1 
47.39 

0.001054  +  7.7941  x  10^29/ 
1.86  + 8.1252  X  10-13/ 
0.9995  -  4.6035  X  10- 1^- 


Signal  2 
35.82 

0.001404-2.8107  X 
2. 147- 4.5929  X 
0.9993  +  3.0037  X  10-11’/ 
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Table  3.3.  Recovery  performance  figures  for  OFDM,  set  2. 


Signal  1 

Signal  2 

Power  Angle  (°) 

45.45 

44.38 

Error  Figure 

0.0004209  -  1.6937  x  lO^^S/ 

0.0004143  + 2.4529  X 

Subspaee  Angle  (°) 

1.176  +  4.1564  X  10-13/ 

1.166- 1.7555  X  10-13/ 

Correlation 

0.9998  -  1.4883  X  10-ii’z 

0.9998  +  6.2367  X  10-12/ 

Table  3.4.  Recovery  performance  figures  for  OFDM,  set  3. 

Signal  1  Signal  2 

Power  Angle  (°)  42.56 

47.23 

Error  Figure 

0.0005003  +  2.2331  x 

0.0004999+1.6457  x 

Subspaee  Angle  (°)  1.282  —  1.9555  x  10  i^/ 

1.281 +3.2425  X  lO-i^/ 

Correlation 

0.9997  +  7.6339  X  10-ii’z 

0.9998  -  1.2653  X  10-12/ 

Table  3.5.  Recovery  performance  figures  for  OFDM,  set  4. 

Signal  1 

Signal  2 

Power  Angle  (°) 
Error  Figure 
Subspaee  Angle  (°) 
Correlation 


45.54 

0.0001- 2.2446  X 

0.573  + 1.1413  xl0-l3/ 

0.9999-1.9919xl0l^z 


42.5 

0.0001001-1.7213x  10-29/ 

0.5732-  1.3908  X 
0.9999  + 2.4284  X  lO-^^i 


Figure  3.1.  OFDM  test  signal,  time  domain. 
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-5  0 
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(a)  First  signal  (b)  Second  signal 

Figure  3.2.  Frequency  domain  of  OFDM  test  transmitted  and  recovered  signals. 
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(a)  Frequency  domain 
Figure  3.3. 


(b)  Time  domain 
Error  of  recovered  OFDM  signals. 
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CHAPTER  4: 

Collection  Geometry  Analysis 


The  recovery  method  in  Section  2.1  depends  upon  an  invertible  A  matrix  in  AS  =  R,  as 
identified  in  Section  2.2.  In  addition  to  an  equal  number  of  transmitters  and  receivers  (m  = 
n)  such  that  A  is  square,  the  collection  geometry  plays  a  significant  role  in  the  invertibility 
of  this  matrix,  as  each  element  of  the  A(/)  block  is  of  the  form  depends 

upon  the  propagation  delay  Tm„.  A  means  to  assess  a  given  collection  geometry  and  an 
initial  analysis  of  the  impact  of  several  such  geometries  upon  the  recovery  performance  is 
provided  in  this  section. 

4.1  Geometry  Assessment 

In  linear  algebra,  the  sensitivity  of  x  to  small  perturbations  in  b  when  solving  the  Ax  =  b 
problem  is  described  by  the  condition  number  of  the  mx  m  A  matrix, 

k{A)  =  —  (4.1) 

where  Oi  is  the  largest  singular  value  and  Om  is  the  smallest  [23].  (Both  are  computed 
via  the  singular  value  decomposition  of  A.)  A  matrix  with  a  small  condition  number  is 
described  as  well- conditioned,  while  a  matrix  with  a  large  condition  number  is  described 
as  ill-conditioned.  A  singular  matrix  is  often  indicated  by  K'(A)  =  oo  [24].  Condition 
numbers  that  yield  acceptable  or  unacceptable  results  are  very  dependent  upon  the  details 
of  the  specific  problem.  See  [23]  or  [24]  for  a  more  complete  discussion  of  the  condition 
number  and  singular  value  decomposition. 

In  the  case  of  (2.7),  the  A  matrix  is  a  block  diagonal  matrix  consisting  of  the  A(z)  blocks 
associated  with  each  frequency  to,  as  defined  in  (2.6).  As  each  of  these  blocks  is  indepen¬ 
dent  of  the  other  blocks  in  the  A  matrix,  it  is  possible  to  assess  the  condition  numbers  of 
the  blocks  independently,  for  the  given  frequency  associated  with  that  block,  rather  than 
simply  evaluating  the  condition  of  the  complete  A  matrix.  This  assessment  provides  in¬ 
sight  as  to  the  recovery  performance  at  each  frequency  used  in  the  recovery  calculation  and 
is  provided  in  conjunction  with  the  frequency-domain  recovery  error  plot  for  each  of  the 
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layouts  presented  in  the  following  section. 


4.2  Layout  Overview  and  Test  Signal 

Four  transmitter  and  receiver  geometries  are  evaluated,  in  addition  to  the  layout  used  in  the 
preceding  sections,  which  is  included  as  a  reference  layout  (Layout  1).  The  transmitter  lo¬ 
cations  remain  constant  throughout  the  analysis.  Layout  2  extends  the  separation  between 
the  receivers  to  500  km  from  the  previous  50  km.  Layout  3  uses  this  same  extended  re¬ 
ceiver  baseline  and  increases  the  altitude  of  the  second  receiver  from  400  km  to  800  km. 
Layout  4  moves  the  second  receiver  such  that  it  appears  at  roughly  five  degrees  above  the 
transmitters’  horizon.  Those  receivers  have  an  altitude  of  200  km.  Lastly,  Layout  5  eval¬ 
uates  the  possibility  of  co-locating  receivers  on  a  single  platform  using  a  1 -meter  receiver 
separation.  These  layout  parameters  are  provided  in  Table  4.1,  with  parameters  that  vary 
from  the  reference  layout  indicated  in  bold  font.  An  illustration  (not  to  scale)  of  the  five 
layouts  is  provided  in  Figure  4.1,  which  includes  a  red  outline  of  the  reference  layout. 

To  facilitate  comparison  across  test  layouts,  a  single  set  of  OFDM  signal  parameters  was 
utilized  for  all  tests.  Those  signal  parameters  (Set  1  from  Section  3.2)  are  provided  in 
Table  4.2.  Additionally,  in  order  to  minimize  any  variation  in  the  recovery  performance 
and  facilitate  comparison  across  the  various  geometries,  the  content  of  each  transmitter’s 
signal  was  maintained  as  a  constant  across  each  of  the  layouts. 

4.3  Results 

Correlation  coefficient  values  for  each  recovered  signal  are  provided  for  the  various  lay¬ 
outs  in  Table  4.3.  The  complete  set  of  metrics  discussed  in  Section  2.3.3  are  provided  in 
Appendix  B.2.  Each  case  has  a  correlation  coefficient  at  or  very  near  one,  and  with  one 
exception,  all  have  imaginary  parts  on  the  order  of  Cmachine-  The  only  exception  is  with 
Layout  5,  with  the  two  receivers  on  the  same  collection  platform,  and  in  that  case  the  imag¬ 
inary  part  is  between  two  and  three  orders  of  magnitude  larger  than  Cmachine  but  remains  on 
the  order  of  10 

Frequency-domain  recovery  error  (as  described  in  Section  2.5.2)  and  condition  number 
plots  are  provided  for  each  layout  in  Figure  4.2  through  Figure  4.6.  The  condition  number 
plots  provide  both  an  overview  of  the  full  spectrum  used  for  recovery,  as  well  as  a  view  of 
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Table  4.1.  Transmitter  and  receiver  geometries,  by  layout.  (In  kilometers.) 


Element 

X 

y 

alt 

Layout  1 

Tx  1 

0 

0 

0 

(Reference) 

Tx2 

0 

30 

0 

Rx  1 

146.2 

0 

400 

Rx2 

146.2 

50 

400 

Layout  2 

Tx  1 

0 

0 

0 

Tx2 

0 

30 

0 

Rx  1 

146.2 

0 

400 

Rx2 

146.2 

500 

400 

Layout  3 

Tx  1 

0 

0 

0 

Tx2 

0 

30 

0 

Rx  1 

146.2 

0 

400 

Rx2 

146.2 

500 

800 

Layout  4 

Tx  1 

0 

0 

0 

Tx2 

0 

30 

0 

Rx  1 

146.2 

0 

400 

Rx2 

1113 

50 

200 

Layout  5 

Tx  1 

0 

0 

0 

Tx2 

0 

30 

0 

Rx  1 

146.2 

50.001 

400 

Rx2 

146.2 

50 

400 

Table  4.2.  Geometry  analysis  test  signal  parameter  summary. 

Value 

Subcarriers  6 

Modulation  Scheme  QPSK 
Bandwidth  (kHz)  600 

Symbols  16 


a  subset  of  the  lower  frequencies,  depicted  on  a  logarithmic  y-axis.  (The  area  contained  in 
the  subset  view  is  indicated  in  the  plot  of  the  full  spectrum  with  a  red  box.)  This  provides 
a  better  understanding  of  the  periodic  nature  of  the  condition  number  that  occurs  with 
some  geometries.  As  the  condition  number  of  the  A(0)  block  is  known  to  be  infinite,  it 
has  been  omitted  from  the  plots.  In  each  layout,  there  is  a  clear  correlation  between  the 
condition  number  at  a  given  frequency  and  the  frequency-domain  recovery  error  depicted 
in  the  corresponding  figure.  It  is  also  worth  noting  that  the  largest  condition  number  (which 
occurs  in  the  Layout  5  case)  is  also  correlated  with  the  largest  recovery  error.  Even  in  that 
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case,  with  a  condition  number  on  the  order  of  10^,  the  frequency-domain  recovery  error 
only  approaches  6  x  and  the  correlation  coefficients  for  that  layout  remain  close  to 

one,  indicating  excellent  recovery  performance. 

The  complete  set  of  signal  and  error  plots  for  each  of  the  layouts  is  provided  in  Ap¬ 
pendix  B.2.  Note  that  Layouts  2,  3,  and  4  involve  larger  propagation  delays  than  Layouts  1 

and  5  and  require  additional  zero  padding  in  the  time  domain.  That  additional  zero  padding 
results  in  samples  at  a  closer  interval  in  the  frequency  domain,  which  is  apparent  in  the 
associated  figures. 

Table  4.3.  Correlation  coefFicient  summary  by  layout  number. 

Signal  1  Signal  2 

Layout  1  0.9987  -  2.6904  x  10^ 0.9991  -  2.496  x  10^ 

Layout  2  0.9994 -b  1.9458  x  0.9995-7.809  x 

Layouts  0.9994-6.9239x  10-1^/  0.9995-5.1021  x  10-11’/ 

Layout4  0.9994-2.9356x  10-1^/  0.9995+ 1.2612  x  lO-i^/ 

Layout  5  0.9987  -  2.2326  x  10  0.9991  +  6.5599  x  IQ-^^i 
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Figure  4.1.  Overview  of  geometry  test  laydowns. 
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Figure  4.2.  Geometry  induced  recovery  errors,  layout  1. 
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Figure  4.3.  Geometry  induced  recovery 


(b)  Condition  number 
errors,  layout  2. 
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(a)  Frequency-domain  error  (b)  Condition  number 

Figure  4.4.  Geometry  induced  recovery  errors,  layout  3. 
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(a)  Frequency-domain  error 

Figure  4.5.  Geometry  induced  recovery 


(b)  Condition  number 
errors,  layout  4. 
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(a)  Frequency-domain  error  (b)  Condition  number 

Figure  4.6.  Geometry  induced  recovery  errors,  layout  5. 
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CHAPTER  5: 
Conclusion 


An  approach  to  signal  recovery  based  upon  the  known  propagation  delays  between  a  set 
of  known  transmitters  and  receivers  was  described  in  this  thesis.  While  results  with  nar¬ 
rowband  signals  do  not  offer  much  initial  promise  for  the  approach,  by  considering  the 
frequency  response  associated  with  delayed  signals  and  the  need  for  the  delay  in  contin¬ 
uous  time  to  resemble  the  circular  shift  associated  with  discrete  time  signals,  admirable 
performance  was  demonstrated  with  the  wideband  signals  typical  of  fourth  generation  cel¬ 
lular  transmissions. 


5.1  Summary  of  Results  and  Identified  Limitations 

Examination  of  recovery  results  with  several  OFDM  signals  as  well  as  in  a  variety  of 
geometries  indicates  that  the  proposed  method  is  worth  further  consideration  and  study. 
The  recovery  approach  produced  correlation  coefficients  approaching  one  for  the  various 
OFDM  signal  parameters  evaluated  as  well  as  the  variety  of  collection  geometries.  The 
use  of  the  condition  number  to  evaluate  prospective  collection  geometries  offers  insight 
into  potential  errors  in  the  frequency  domain  of  the  recovered  signals  and  the  viability  of  a 
potential  collection  from  a  geometry  perspective. 

A  number  of  limitations  exist  that  must  be  addressed  before  the  approach  is  viable  for  real- 
world  signal  recovery.  Performance  in  the  presence  of  noise  is  a  considerable  limitation. 
The  above  results  indicate  the  need  for  an  additional  35  dB  of  SNR,  a  considerable  margin 
under  any  realistic  circumstances.  The  known  singularity  at  /  =  0  Hz  may  also  be  a  limi¬ 
tation  if  applied  to  a  signal  that  has  a  non-zero  DC  subcarrier.  Computational  complexity 
must  also  be  considered,  given  the  potential  size  of  the  A  matrix.  These  limitations  and 
potential  solutions  are  further  described  in  the  following  section. 

5.2  Future  Work 

There  are  a  number  of  areas  that  merit  additional  study,  either  to  address  identified  short¬ 
comings  of  the  approach  as  it  currently  exists  or  to  extend  the  work  in  order  to  ensure 
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adequate  performance  in  the  intended  application. 


5.2.1  Received  Signal  Model  Improvement 

As  discussed  in  Section  2.5.1,  the  model  currently  constructs  wideband  signals  as  the  sum 
of  delayed  copies  of  time-limited  transmissions.  For  small  differences  in  delay  times,  this 
results  in  a  small  fraction  of  each  received  signal  that  only  contains  signal  from  a  single 
transmitter,  as,  due  to  differences  in  the  propagation  delays  between  the  transmitter  and 
receiver  pairs,  the  signals  from  the  other  transmitters  have  either  not  yet  arrived  or  have 
completely  arrived.  For  longer  propagation  delays  coupled  with  shorter  signals,  this  can 
potentially  result  in  the  signals  arriving  at  the  receiver  at  distinctly  different  times. 

While  the  recovery  process  does  not  take  advantage  of  this  difference  in  arrival  times,  it 
is  more  likely  that  a  base-station  transmitter  would  exhibit  a  much  larger  transmit  duty 
cycle  and  a  more  accurate  model  would  ensure  that  the  each  sample  of  the  received  signal 
contained  a  contribution  from  each  of  the  transmitters.  This  would  require  the  generation 
of  longer  duration  signals  and  the  application  of  a  windowing  or  indicator  function  such 
as  (2.26),  where  the  indicator  domain  A  is  limited  to  a  period  of  time  when  all  signals  are 
present  in  the  receiver. 

5.2.2  Computational  Complexity 

In  (2.7),  the  A  matrix  is  niN  x  nN,  where  m  and  n  are  the  numbers  of  receivers  and  trans¬ 
mitters,  respectively.  The  parameter  N  is  the  number  of  discrete  frequencies,  which  is  a 
function  of  the  desired  signal  duration,  sampling  rate  and  amount  of  zero  padding  required 
for  the  given  signal  delay.  Increasing  the  sample  rate  in  order  to  accommodate  higher 
bandwidth  signals  or  the  extending  the  buffer  length  when  zero  padding  signals  with  larger 
propagation  delays  has  the  potential  to  dramatically  increase  the  size  of  the  A  matrix  and 
the  computation  time  associated  with  signal  recovery  as  a  result.  Of  note,  in  this  work, 
this  limitation  required  a  lower  altitude  for  the  off-nadir  receiver  in  Layout  4  of  the  geom¬ 
etry  test.  It  may  be  possible  to  use  the  difference  in  propagation  time  across  the  various 
transmitter/receiver  pairs  vice  the  absolute  time,  which,  in  some  geometries  would  offer  a 
considerable  savings  in  the  amount  of  zero  padding  required. 

Some  improvement  may  also  be  available  as  a  result  of  the  fact  that  the  A  matrix  is  par¬ 
ticularly  sparse.  Use  of  MATLAB’s  sparse  matrix  tools  should  both  decrease  memory 
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requirements  and  decrease  computation  time  given  that  the  number  of  non-zero  entries  in 
the  A  matrix  grows  as  0{N)  while  size  grows  as  0{N^).  Similarly,  as  currently  defined,  it 
is  possible  to  break  the  AS  =  R  equation  into  a  set  of  A(/)S(/)  =  R(/)  equations  for  each 
discrete  frequency  i  and  solve  those  independently.  This  would  complicate  approaches  that 
leverage  any  potential  correlation  between  adjacent  frequencies  to  improve  performance 
with  noise  or  in  the  m  <  n  case.  (See  below.) 

5.2.3  Accuracy  of  Delay  Knowledge 

The  preceding  work  is  entirely  dependent  upon  and  assumes  completely  accurate  knowl¬ 
edge  of  the  delays  associated  with  signal  propagation  between  each  transmitter  and  receiver 
pair  and  uses  double  precision  variables  for  the  associated  computations.  There  is  no  error 
included  in  any  transmitter  or  receiver  position  data,  either  when  constructing  the  received 
signals  or  when  solving  for  the  signals  as  originally  transmitted. 

Such  levels  of  accuracy  are  not  available  when  using  real-world  transmitters  and  receivers, 
and  it  is  necessary  to  evaluate  performance  under  varying  accuracy  conditions.  That  analy¬ 
sis  should  likely  include  both  an  analytical  evaluation  of  the  pertinent  recovery  equations  to 
quantify  the  impact  of  a  small  change  in  the  expected  signal  delay  on  the  recovered  signals 
as  well  as  a  Monte  Carlo  simulation  with  varying  error  levels  to  confirm  that  analysis. 

The  matter  is  complicated  by  the  interactions  between  the  different  delays  in  the  A  matrix. 
If  each  transmitter  and  receiver  has  an  associated  position  error,  the  geometric  arrangement 
of  those  antennas  will  have  a  different  effect  on  the  delay  between  each  transmitter  and 
receiver.  Those  various  errors  are  then  combined  in  the  A  matrix  and  impact  the  recovered 
signal.  It  would  be  useful  to  understand  the  impact  of  a  positional  error  associated  with  a 
single  transmitter  or  receiver  as  well  as  a  set  of  comparable  errors  across  the  complete  set 
of  transmitters  and  receivers. 

5.2.4  Noise 

As  identified  in  Section  2.6,  noise  poses  a  considerable  problem  for  the  recovery  process  as 
it  solves  the  AS  =  R  equation  while  AS  -1-  r]  =  R  is  a  more  accurate  model  of  the  received 
signals.  One  potential  approach  to  mitigating  the  impact  of  the  noise  on  signal  recovery 
relies  upon  the  fact  that  the  signals  of  interest  are  transmitting  digital  data  in  accordance 
with  published  standards,  such  as  [7]. 
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Those  standards  include  a  specific  symbol  library,  and  forcing  the  recovery  signal  to  ad¬ 
here  to  the  library  may  mitigate  some  of  the  impacts  of  noise.  After  an  initial  attempt 
at  signal  recovery,  the  recovered  signals  might  be  processed  as  they  would  by  a  typical 
commercial  receiver,  including  the  detection  process,  yielding  channel  symbols  [25].  By 
re-encoding  those  recovered  channel  symbols,  the  resulting  signal  is  effectively  a  clean 
version  of  the  recovered  signal  (precisely  what  was  transmitted,  in  the  ideal  case),  and  any 
noise-induced  differences  between  the  recovered  symbol  and  the  re-encoded  symbol  are 
eliminated.  While  this  does  not  protect  against  any  difference  that  effectively  results  in  an 
incorrect  channel  symbol,  it  may  offer  some  performance  improvement  over  the  existing 
recovery  process.  While  it  likely  involves  a  considerable  increase  in  complexity,  it  may  be 
possible  to  extend  this  decoding/re-encoding  process  to  leverage  any  error  correction  cod¬ 
ing  present  in  the  signal-of-interest,  further  mitigating  the  impact  of  noise  on  the  recovered 
signal. 

It  is  also  important  to  recognize  that  the  condition  number  serves  as  an  indicator  of  the 
sensitivity  of  the  recovered  signal  to  small  perturbations  in  the  received  signals.  When 
the  condition  number  is  larger,  the  noise  in  the  received  signals  will  have  a  greater  impact 
on  the  recovered  signal  than  in  those  cases  with  smaller  condition  numbers.  With  further 
study,  it  may  be  possible  to  develop  an  empirical  understanding  of  the  cases  when  a  larger 
condition  number  requires  some  alternative  recovery  approach  or  noise  mitigation  strategy 
prior  to  or  as  part  of  the  recovery  process. 

5.2.5  Evaluation  of  Recovered  Signals  as  Bistatic  Illuminators 

Ultimately,  the  signal  recovery  process  is  intended  to  support  bistatic  SAR  imaging,  with 
a  recovered  signal  serving  as  the  illuminator  for  the  imaging  process.  It  is  important  to 
understand  the  implications  to  the  imaging  process  when  a  recovered  signal  is  utilized  vice 
the  more  typical  case  employing  a  known  pulse  shape.  Holder  [18]  developed  a  model  of  a 
monostatic  SAR  system  which  could  be  extended  to  accommodate  bistatic  signals  and  then 
used  to  evaluate  the  impact  of  the  recovery  process. 

There  are  two  fundamental  changes  to  the  existing  model  necessary  to  support  that  evalua¬ 
tion.  First,  the  model  needs  to  be  extended  to  accommodate  bistatic  imaging,  including  the 
incorporation  of  bistatic  delay  calculations  for  the  layout  of  targets.  In  the  monostatic  case. 
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the  round-trip  fast-time  delay  is  [2] 


2^/x^  +  {y  -  u)^ 
c 


(5.1) 


where  x  and  y  deseribe  the  target  position,  u  is  the  radar’s  position  along  the  y-axis,  and  c  is 
the  speed  of  light.  In  the  bistatie  case,  this  delay  calculation  must  account  for  the  location 
of  the  transmitter  and  is  no  longer  a  round  trip  time.  In  that  case,  the  delay  is  given  by 


t 


-  \l {x- pf  +  {y - qY  +  \l +  {y - uf 


(5.2) 


where  p  and  q  describe  the  location  of  the  transmitter.  These  delay  calculations  are  then 
used  to  construct  the  signals  received  during  the  collection  process. 


The  SAR  model  described  by  Holder  uses  Spatial  Frequency  Interpolation  for  the  image 
reconstruction  process,  which  converts  the  received  signals  to  a  recognizable  image.  This 
algorithm  requires  two-dimensional,  frequency-domain  matched  filtering  [2] 


S{(0,ku)Sl{(0,ku),  (5.3) 

where  5(tt),A:„)  is  the  frequency-domain  representation  of  the  received  SAR  signal  and 
SQ{(0^ku)  is  the  two-dimensional  matched  filter  frequency  response.  While  Soumekh  [2] 
describes  the  derivation  of  SQ{G),ku)  for  the  monostatic  case,  it  will  be  necessary  to  adapt 
that  to  the  bistatic  case  in  order  to  create  the  SAR  images  required  to  properly  evaluate 
signal  recovery  performance.  The  recovered  signal  is  used  directly  in  that  matched  filtering 
process  and  is  referred  to  as  the  reference  signal  in  both  [2]  and  [18]. 

The  ability  to  study  the  impact  of  using  recovered  signals  for  SAR  image  formation  then  al¬ 
lows  further  examination  of  the  impact  of  noise  and  positional  errors  on  the  overall  system. 
While  both  of  those  can  be  worked  somewhat  independently,  that  work  ultimately  impacts 
the  quality  of  the  SAR  image,  and  understanding  the  impacts  of  both  noise  and  positional 
error  on  the  image  quality  is  critical  to  the  utility  of  the  approach. 
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5.2.6  Adaptation  to  m<n  Scenarios 

Given  the  proliferation  of  LTE  and  WiMax  transmitters,  eoupled  with  the  fiseal  eonstraints 
and  limited  resouree  availability,  a  signal  reeovery  approach  that  functions  when  there  are 
more  transmitters  than  receivers  is  particularly  desirable.  This  is  a  classical  example  of  an 
underdetermined  problem,  as  given  the  resulting  shape  of  the  A  matrix,  there  is  not  enough 
information  available  to  yield  a  single  solution  to  the  AS  =  R  matrix  equation.  Common 
approaches  to  this  underdetermined  problem  such  as  seeking  the  solution  that  contains  the 
most  zero  elements  or  that  minimizes  an  error  function  (often  referred  to  as  the  method  of 
Least  Squares)  do  not  seem  well  suited  to  this  particular  problem,  as  there  is  no  expectation 
that  the  original  waveforms  adhered  to  either  of  those  paradigms. 

By  oversampling  the  received  signals,  it  is  possible  to  develop  some  redundant  informa¬ 
tion  in  those  signals.  In  the  frequency  domain,  this  means  that  the  adjacent  samples  are  no 
longer  independent.  By  incorporating  that  assumption  into  the  construction  of  the  A  ma¬ 
trix,  such  that  the  solution  for  a  given  point  in  the  frequency-domain  representation  of  the 
signal  is  dependent  upon  other  surrounding  values,  it  may  be  possible  to  develop  a  better 
approximation  for  the  original  signal. 


50 


APPENDIX  A: 
Signal  Interpolation 


In  order  to  evaluate  the  signal  isolation  approach  described  in  Section  2.1,  it  is  necessary  to 
generate  a  collection  of  signals  that  contains  a  set  of  original  signals  at  various  delays.  This 
is  straightforward  in  the  narrowband  case  which  uses  signals  defined  over  a  time  limited 
period  t  of  the  form 

Snit)  =  AnCO&{a>nt)  (A.l) 

as  the  cosine  is  a  continuous  function  and  the  required  delays  can  be  computed  as 
Sn{t  —  Xmn)-  In  the  wideband  case,  which  simply  uses  time-limited  AWGN  as  the  signal-of- 
interest,  one  approach  to  computing  the  various  delays  is  to  interpolate  the  original  signal 
at  the  necessary  delayed  sample  times.  This  interpolation  is  commonly  performed  using 
Ideal  Low  Pass  Filter  Interpolation.  As  an  alternative,  an  approach  utilizing  a  phase  shift 
in  the  frequency  domain  was  evaluated.  Both  are  described  below  with  a  comparison  of 
results  of  the  two  approaches. 


A.l  Ideal  Low  Pass  Filter  Interpolation 

Interpolation  allows  recovery  of  information  from  the  original  continuous  time  signal  that 
is  not  immediately  available  after  the  sampling  process.  This  approach  to  interpolation 
utilizes  a  low  pass  filter  to  counter  the  effects  of  the  sampling  process.  The  sampling 
process  and  its  effects  in  the  time  and  frequency  domain  are  reviewed  below  prior  to  an 
explanation  of  interpolation  using  an  ideal  low  pass  filter. 


A.  1.1  Sampling  in  the  Time  Domain 

Ideal  sampling  of  a  continuous  time  signal  s{t)  can  be  considered  as  the  product  of  that 
signal  with  a  train  of  unit  impulses 

-j-oo 

P{t)=  d{t-nT),  (A.l) 
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yielding  the  sampled  signal 


=x{t)p{t) 

-\-oo 

=  ^  x{t)5{t  —  tiT) 

n=—co 

-\-oo 

=  Y,  x{nT)d(t-nT)  (A.3) 

n=—oo 

as  a  train  of  pulses  with  area  determined  by  the  signal  at  each  sample  interval.  For  con¬ 
venience,  the  notation  x[n\  is  used  to  represent  a  discrete  time  equivalent  of  x{nT).  This 
process  is  illustrated  in  Figure  A. la.  Once  sampled,  the  continuous  time  information  is  no 
longer  directly  available  —  interpolation  is  required  to  recover  signal  information  at  other 
than  the  sampled  instants. 


Xs{t) 


1 


± 


Uf) 

A/T 


-2T  -TO  T  2T  3T  4T  5T  -2/.  -fs  -fi  0  /i  /, 
(a)  Time  domain  (b)  Frequency  domain 

Figure  A.l.  Sampling  a  continuous  time  signal,  after  [19]. 


2A 
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A.1.2  Frequency-domain  Perspective 

In  the  frequency  domain,  depicted  in  Figure  A.  lb,  the  impulse  train  used  for  sampling  is 
also  an  impulse  train,  and  taking  /^  =  1/r,  is  described  by 

-)-CX3 

P{f)=fs  L  5{f-kf,).  (A.4) 

k=  —  aa 

The  multiplication  of  s{t)  and  p{t)  in  the  time  domain  corresponds  to  convolution  in  the 
frequency  domain,  such  that 

Uf)=x{f)*Pif) 

-\-oo 

=fs  L  x{f)*d{f-kf,) 

k=  —  ca 
-\-oo 

=  /.  L  X{f-kf,),  (A.5) 

k=—oo 

which  is  a  train  of  the  original  spectrum  repeated  at  multiples  of  fs  in  the  frequency  domain 
[19],  [21].  As  long  as  the  sampling  rate  exceeds  that  required  by  the  Nyquist  sampling 
theorem,  aliasing  is  avoided. 


Mf) 


-2fs  -fs-h  0  /i  /,  If, 


Figure  A. 2.  Ideal  low  pass  filter. 


A.1.3  Interpolation  in  the  Frequency  Domain 

In  the  absence  of  aliasing,  in  order  to  recreate  the  original  continuous  time  signal,  the 
original  frequency  spectrum  can  be  obtained  by  applying  an  ideal  low  pass  filter  to  isolate 
the  spectrum  of  the  original  signal,  as  depicted  in  Figure  A. 2  [12].  The  ideal  low  pass  filter 
is  described  by 

Hilpf(/)  =  (4)  (A.6) 
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where 


n( 


jc  = 


1  if|x|<l/2 

0  otherwise 

is  known  as  the  reetangle  function.  Applying  the  ideal  low  pass  filter  to  (A. 5)  gives 


(A.7) 


y(f)  =  ^ilpf(/)A,(/)  =  Kf,X(f),  (A.8) 

and  selecting  K  =  1/fs  gives  Y  (/)  =X{f).  The  inverse  Fourier  transform  can  then  be  used 
to  recover  the  original  time-domain  signal. 


A. 1.4  Interpolation  in  the  Time  Domain 


Figure  A. 3.  Signal  interpolation  as  a  sum  of  scaled  and  shifted  sinc()  pulses. 

The  multiplication  of  Xs{f)  by  the  ideal  low  pass  filter’s  frequency  response  in  (A.8)  is 
convolution  in  the  time  domain,  such  that  y{t)  =  /?iLFp(t)  Substituting  (A. 3)  then 

gives 


n 

=  ^A:[n]/ziLFp(0  *  5  (t-nT) 

n 

=  ^.r[n]hiLFp(f-«T).  (A.9) 

n 
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Given  the  Fourier  transform  relationship 


(A.IO) 

and  ehoosing  B  =  fs/2  and  K  =  l/fs,  (A.9)  beeomes 

^(0  =  (A.ll) 

n 

whieh  is  a  sum  of  shifted  sine()  pulses,  each  scaled  by  a  corresponding  coefficient  from 
the  x[n]  vector  [12].  This  is  illustrated  in  Figure  A. 3,  which  depicts  the  constituent  sinc() 
functions  that  are  added  to  construct  the  x{t)  signal.  Each  sinc()  pulse  is  zero  at  all  mul¬ 
tiples  of  T  other  than  the  peak,  and  each  sample  point  only  receives  a  contribution  from  a 
the  single  sinc()  that  corresponds  to  that  particular  x[n]  value.  (Note  that  all  but  one  of  the 
sinc()  curves  are  zero  below  each  of  the  red  circles  corresponding  to  a  sample  in  x[n].)  All 
other  points  in  the  x{t)  signal  are  composed  of  contributions  from  each  of  the  sinc()  pulses. 
MATLAB  code  to  perform  this  interpolation  function  is  provided  in  Appendix  C.6  as  the 
function  interp_time  .m.  As  MATLAB  works  with  sampled  signals,  this  code  returns  the 
signal  resampled  at  new  time  values,  which  are  specified  as  an  argument  to  the  function. 


A.1.5  Error  Analysis 

As  a  baseline  for  later  sections,  it  is  useful  to  understand  the  amount  of  error  induced  by 
this  interpolation  process.  As  a  challenging  case,  a  signal  x{t)  consisting  of  AWGN  was 
sampled  to  create  x[n].  A  sampled  version  of  x{t  —  t)  was  generated  by  interpolating  s[n] 
at  a  delayed  point  in  time,  where  x  is  an  integer  multiple  of  the  sampling  interval  and  a 
significant  fraction  of  the  overall  signal  duration  to  make  the  delay  apparent  in  the  follow¬ 
ing  figures.  For  comparison,  a  circular  shift  corresponding  to  x  was  applied  to  the  same 
sampled  AWGN  signal.  Those  signals  are  presented  in  Figure  A.4,  along  with  the  differ¬ 
ence  between  those  two  signals  (plotted  on  a  separate  y-axis  to  emphasize  the  dramatic 
difference  in  scale.)  The  total  energy  in  both  shifted  signals  was  computed  as 

Ei  =  'S.{  XiTsamp?  (A. 12) 
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where  Xi  is  a  eolumn  veetor  of  the  signal-of-interest  with  sampling  interval  Tsamp,  and  Xi* 
is  the  eonjugate  transpose  of  the  column,  such  that  Xi*Xi  computes  the  inner  product  of  the 
signal.  This  was  compared  to  the  energy  in  the  original  x[n].  The  differences  in  energy  are 
presented  in  Table  A.  1  and  show  error  on  the  order  of  £machine>  which  indicates  that  shifting 
via  interpolation  is  a  viable  option. 


0  500  1000  1500  2000  2500  3000  3500  4000  4500 

Time  sample 

Figure  A. 4.  Signals  and  error  from  shifting  via  interpolation. 


Table  A.l.  Energy  summary  for  shifting  via  interpolation. 

Energy  Error  (J) 

Circular  Shift  4.3368  x  lO-i*^ 
Interpolation  2.4286  x  10^^^ 


A.2  Fourier  Phase  Shifting  for  Interpolation 

In  the  continuous  time  domain,  the  Eourier  transform  has  a  property  that  relates  a  shift  in 
the  time  domain  to  a  phase  shift  in  the  frequency  domain  [19],  such  that 

x{t)^X{(0)  (A.13) 

x{t-td)^X{o})e-^^^‘^.  (A.14) 

In  the  discrete  time  domain,  the  discrete  Eourier  transform  has  a  similar  property  such 
that  a  circular  shift  in  the  discrete  time  domain  results  in  a  phase  shift  in  the  frequency 
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domain  [21].  That  is, 


(A.15) 


(A.16) 


where  n  is  the  index  or  sample  number,  ©  represents  the  left  eireular  shift  operation,  m 
is  the  integer  number  of  samples  to  shift  by,  N  is  the  total  number  of  samples,  and  k  is 
the  discrete  frequency  index.  With  adequate  zero  padding,  this  circular  shift  resembles  the 
time  delay  from  the  continuous  time  domain,  where  the  delay  is  given  by  =  mTsamp  and 
Tsamp  is  the  sampling  interval.  This  allows  for  shifting  of  signals  by  integer  multiples  of  the 
sampling  interval  by  applying  the  FFT,  multiplying  by  the  appropriate  phase  shift  in  the 
frequency  domain,  and  returning  to  the  time  domain  with  the  IFFT.  While  this  convenient 
for  signals  that  are  already  in  the  frequency  domain,  it  is  likely  excessive  effort  for  a  signal 
in  the  time  domain,  as  a  circular  shift  operation  is  accomplished  with  a  straightforward 
manipulation  of  the  index  numbers. 

For  comparison  to  other  cases,  plots  of  the  delayed  signals  and  associated  error  between 
this  method  and  the  circular  shift  in  the  time  domain  are  provided  in  Figure  A. 5.  As  this 
method  involves  multiplication  by  complex  terms,  plots  of  both  the  real  and  imaginary  parts 
of  the  signal  are  provided.  Similarly,  the  energy  difference  between  the  phase  shifted  and 
circularly  shifted  signals  are  provided  in  Table  A. 2,  which  also  includes  entries  for  both  the 
real  and  imaginary  portions  of  the  error.  Each  is  well  below  Cmachine- 


Real  part  of  time  domain  signals 


-2 


-4 


2 


CircshiftO 
Phase  Shift 
Difference 


0  E 


-2 


500  1000  1500  2000  2500  3000  3500  4000 

Time  sample 

Imag  part  of  time  domain  signals 


4500 


10' 


_  0.5' 


-0.5' 


-2 


-1 


— ^ ^ ^ ^ ^ ^ ^ ^ - ^_4 

500  1000  1500  2000  2500  3000  3500  4000  4500 

Time  sample 


Figure  A. 5.  Signals  and  error  from  integer  phase  shifts  in  the  frequency  domain. 
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Table  A. 2.  Energy  summary  for  integer  phase  shifts  in  the  frequency  domain. 

Energy  (J) 

Error  8.9355  x 

Real(Error)  8.934  x  10^29 
Imag(Error)  E51  x  10^^^ 

A.2.1  Fractional  Delays 

The  opportunity  to  use  non-integer  or  fraetional  delays  for  shifting  is  appealing  as  it  ae- 
eomplishes  an  interpolation  to  get  the  new  set  of  values  of  x{t)  that  are  required  but  not 
available  direetly  in  x[n] .  The  plots  of  the  signals  and  error  that  result  from  sueh  an  attempt 
are  provided  in  Eigure  A. 6.  As  the  eireular  shift  does  not  aeeommodate  fraetional  shift 
amounts,  an  interpolation  was  used  as  the  baseline  signal  for  eomparison  and  when  eom- 
puting  both  the  error  plot  and  energy  errors,  whieh  are  presented  in  Table  A. 3.  The  error 
levels  in  this  ease  are  mueh  more  signifieant  than  those  with  integer  shift  amounts,  and, 
more  signifieantly,  the  error  in  the  imaginary  portion  of  the  signal  is  an  order  of  magnitude 
larger  than  that  of  the  real  portion.  Intuitively,  as  the  proeess  begins  with  a  real  signal  and 
the  intent  is  to  merely  delay  the  signal  in  time,  the  result  should  not  have  an  imaginary 
portion. 


Real  part  of  time  domain  signais 


Time  sample 

Figure  A. 6.  Signals  and  error  from  fractional  phase  shifts  in  the  frequency  domain. 

Table  A. 3.  Energy  summary  for  fractional  phase  shifts  in  the  frequency  domain. 

Energy  (J) 

Error  6.1483  x  10^'' 

Real(Error)  7.2528  x  10^^ 

Imag(Error)  5.423  x  10“^ 
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Error  Source 

The  Fourier  transform,  in  both  the  continuous  and  discrete  forms,  has  a  number  of  prop¬ 
erties  that  are  helpful  when  considering  this  error.  In  particular,  when  x{t)  is  real,  the 
frequency-domain  description  of  the  same  signal  X{f)  will  be  conjugate  symmetric,  such 
that  X{—f)  —  X*{f)  [19].  The  inverse  is  also  true,  such  that  a  signal  must  be  conjugate 
symmetric  in  the  frequency  domain  in  order  to  be  real  in  the  time  domain. 

Recall  that  in  the  case  of  discrete  signals,  the  frequency  domain  is  periodic.  In  that  case, 
conjugate  symmetry  requires  that  X[n]  =  X*[N  —  n]  where  the  signal  is  periodic  over  N 
points,  n  <  N,  and  the  first  point  is  n  =  0.  (In  MATLAB,  indices  begin  at  n  =  I,  so  in  that 
case  conjugate  symmetry  requires  X[n]  =  X*[N  +  2  —  n].  The  remainder  of  this  discussion 
will  assume  the  first  index  is  n  =  0.)  Two  of  these  points  are  significant  —  n  =  0  and 
n  =  N /2.  When  n  =  0,  the  above  relationship  requires  conjugate  symmetry  with  the  point 
at  n=N,  which  exceeds  the  length  of  the  sampled  signal  but  corresponds  to  itself  when  the 
periodic  nature  of  the  DFT  frequency  domain  is  considered.  The  same  is  true  for  n  =  N/2, 
which  must  also  be  conjugate  symmetric  with  itself,  or,  more  simply,  real. 

When  shifting  a  signal  by  a  fraction  of  the  sampling  interval,  the  term  must  be 

conjugate  symmetric  in  order  for  the  resulting  signal  to  be  real-valued  in  the  time  domain. 
That  requires  that  both  the  n  =  0  and  n  =  N/2  values  in  that  term  to  be  real  valued,  which 
is  not  the  case  when  shifting  by  fractions  of  a  sampling  interval.  In  that  case  the  n=N /2 
value  will  be  complex,  and  the  resulting  signal  will  have  an  imaginary  part  in  the  time 
domain. 


Forcing  Conjugate  Symmetry 

As  a  real  valued  result  is  expected,  the  n  =  N /2  point  can  be  artificially  forced  to  a  real 
value  by  eliminating  the  imaginary  portion  of  that  sample  prior  to  applying  the  term. 
In  that  case,  the  phase  shift  term  is  conjugate  symmetric,  and  the  resulting  signal  should  be 
shifted  in  the  time  domain  without  incurring  an  imaginary  portion.  Results  of  that  process 
are  provided  in  Figure  A.7,  which  again  depicts  errors  relative  to  a  signal  delayed  via 
interpolation.  The  error  in  the  signal  energy  is  presented  in  Table  A.4,  which  shows  the 
energy  in  the  imaginary  portion  has  been  reduced  to  zero  as  expected. 
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Real  part  of  time  domain  signais 


Figure  A. 7.  Signals  and  error  from  fractional  phase  shifts  with  forced  symmetry. 

Table  A. 4.  Energy  summary  for  fractional  phase  shifts  with  forced  symmetry. 

Energy  (J) 

Error  7.2528  x  10^^ 

Real(Error)  7.2528  x  10^^ 

Imag(Error)  0 

A.3  Computation  Time  Comparison 

If  speed  is  a  eonsideration,  and  the  applieation  ean  tolerate  the  antieipated  error,  delay 
by  the  fraetional  phase  shift  method,  ineluding  a  correction  to  eliminate  the  imaginary 
portion  of  the  result,  should  be  considered.  The  interpolation  method  requires  repeated 
computation  of  the  sinc()  function,  which,  for  a  large  number  of  samples  in  the  signal  or 
desired  points  of  interpolation,  can  be  computationally  intensive.  The  phase  shift  method 
requires  computation  of  the  EET  and  lEET.  A  comparison  of  the  times  required  for  the 
above  evaluation  is  provided  in  Table  A. 5.  While  this  is  just  a  small  data  set  for  a  single 
signal,  it  is  clear  that  the  interpolation  requires  considerably  more  time  to  compute  than  the 
delays  via  phase  shift  in  the  frequency  domain.  Eittle  more  should  be  extracted  from  this 
data,  given  the  minimal  number  of  runs  to  evaluate  performance. 


A.4  Summary 

When  shifting  an  integer  number  of  time-domain  samples,  a  circular  shift  operation  (e.g. 
using  circshiftO  in  MATEAB)  is  a  straightforward  option,  as  this  simply  involves  ad¬ 
justing  the  index  numbers  of  the  sampled  data  and  should  yield  minimal,  if  not  zero,  error. 
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Table  A. 5.  Calculation  times  for  various  shifting  methods. 

Calculation  time  (ms) 


Circular  Shift 

3.3 

Interpolation  (integer) 

355.2 

Interpolation  (fractional) 

358.2 

Phase  Shift  (integer) 

1.7 

Phase  Shift  (fractional) 

2.0 

Phase  Shift  (fractional  w/  correction) 

4.4 

For  fractional  shifts,  options  include  either  interpolation  at  the  delayed  time  or  use  of  a 
fractional  phase  shift  in  the  frequency  domain.  As  discussed  above,  there  is  some  error 
between  those  two  methods,  and  while  a  true  “control”  for  that  error  is  not  available,  it  is 
reasonable  to  attribute  the  discrepancy  to  the  fractional  phase  shift  method,  as  it  involves 
a  departure  from  the  typical  case  while  interpolation  accuracy  should  not  be  affected  by 
the  times  at  which  the  signal  is  interpolated.  Depending  upon  the  application,  computation 
time  might  be  a  considerable  argument  for  use  of  the  delay  via  phase  shift  in  the  frequency 
domain. 
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APPENDIX  B: 

Additional  Analysis  Tables  and  Figures 


B.l  OFDM  Signal  Recovery 
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Figure  B.l.  OFDM  test  signal,  time  domain,  parameter  set  1. 


Frequency  (kHz) 


Frequency  (kHz) 


(a)  First  signal  (b)  Second  signal 

Figure  B.2.  Frequency  domain  of  OFDM  test  transmitted  and  recovered  signals,  parameter  set 
1. 
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(a)  Frequency  domain  (b)  Time  domain 

Figure  B.3.  Error  of  recovered  OFDM  signals,  parameter  set  1. 
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Figure  B.4.  OFDM  test  signal,  time  domain,  parameter  set  2. 
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(a)  First  signal 

Figure  B.5.  Frequency  domain  of 
2. 
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(b)  Second  signal 

OFDM  test  transmitted  and  recovered  signals,  parameter  set 


(a)  Frequency  domain 
Figure  B.6.  Error  of 


(b)  Time  domain 
recovered  OFDM  signals,  parameter  set  2. 
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Figure  B.7.  OFDM  test  signal,  time  domain,  parameter  set  3. 
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Figure  B.8. 
3. 


Frequency  domain  of  OFDM  test  transmitted  and  recovered  signals. 


parameter  set 
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(b)  Time  domain 


Figure  B.9.  Error  of  recovered  OFDM  signals,  parameter  set  3. 
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Figure  B.IO.  OFDM  test  signal,  time  domain,  parameter  set  4. 
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Frequency  (kHz) 


Frequency  (kHz) 


Frequency  (kHz)  Frequency  (kHz) 


(a)  First  signal  (b)  Second  signal 

Figure  B.ll.  Frequency  domain  of  OFDM  test  transmitted  and  recovered  signals,  parameter  set 
4. 


(a)  Frequency  domain  (b)  Time  domain 

Figure  B.12.  Error  of  recovered  OFDM  signals,  parameter  set  4. 
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B.2  Geometry  Analysis 


Table  B.l.  Recovery  performance  figures  for  layout  1. 


Signal  1 

Signal  2 

Power  Angle  (°) 

39.17 

48.65 

Error  Figure 

0.002572-9.5719x  10-29/ 

0.001851+9.2041  x  lO-^^/ 

Subspace  Angle  (°) 

2.907  +  3.0395x  10-14/ 

2.466  + 3.3244  X  10-13/ 

Correlation 

0.9987  -  2.6904  X  10-12/ 

0.9991- 2.496  X  10  i^/ 

Table  B.2.  Recovery  performance  figi 

ures  for  layout  2. 

Signal  1 

Signal  2 

Power  Angle  (°) 

39.11 

48.65 

Error  Figure 

0.001286- 3.2622  X 

0.0009253  +  1.4881  x 

Subspace  Angle  (°) 

2.055-3.1089x  10-13/ 

L743  + 1.4709  X 

Correlation 

0.9994+ 1.9458  X  10 

0.9995  -  7.809  X  10-10/ 

Table  B.3.  Recovery  performance  figi 

ures  for  layout  3. 

Signal  1 

Signal  2 

Power  Angle  (°) 

39.11 

48.65 

Error  Figure 

0.001286- 2.6708  X 

0.0009253  +  6.2858  X  10-28/ 

Subspace  Angle  (°) 

2.055  +  1.1063x  10-13/ 

1.743  +  9.6102  X  10-13/ 

Correlation 

0.9994  -  6.9239  X  10-12/ 

0.9995-5.1021  x 

Table  B.4.  Recovery  performance  figi 

ures  for  layout  4. 

Signal  1 

Signal  2 

Power  Angle  (°) 

39.11 

48.65 

Error  Figure 

0.001286-2.9969x  10-30/ 

0.0009253  +  2.2486  x  10-29/ 

Subspace  Angle  (°) 

2.055  +  4.6904x  10-14/ 

1.743- 2.3755  X  10-13/ 

Correlation 

0.9994  -  2.9356  X  10-12/ 

0.9995  +  1.2612  X  10-10/ 
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Table  B.5.  Recovery  performance  figures  for  layout  5. 


Power  Angle  (°) 
Error  Figure 
Subspaee  Angle  (°) 
Correlation 


Signal  1 
39.17 

0.002572  +  5.2147  X  10-23/ 

2.907  + 2.5224  X  10-‘^/ 

0.9987  -  2.2326  X  10-12/ 


Signal  2 
48.65 

0.001851  -3.0583  X  10-21/ 
2.466-8.7371  x  10-i°/ 
0.9991  +  6.5599  X  10  13/ 


«'«  ’f=i’  /a  r  .'fe 


(a)  Frequency-domain  error 

Figure  B.13.  Geometry  induced 


(b)  Condition  number 
recovery  errors,  layout  1. 
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(b)  Condition  number 


Figure  B.14.  Geometry  induced  recovery  errors,  layout  2. 
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(a)  Frequency-domain  error  (b)  Condition  number 

Figure  B.15.  Geometry  induced  recovery  errors,  layout  3. 


(a)  Frequency-domain  error 

Figure  B.16.  Geometry  induced 


(b)  Condition  number 
recovery  errors,  layout  4. 
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(a)  Frequency-domain  error 

Figure  B.17.  Geometry  induced 


(b)  Condition  number 
recovery  errors,  layout  5. 


71 


ZS,  (deg)  ,g  I  ZS,  (deg) 


Frequency  (kHz) 


Frequency  (kHz) 


(a)  First  signal  (b)  Second  signal 

Figure  B.18.  Frequency  domain  of  recovered  signals,  layout  1. 


(a)  First  signal  (b)  Second  signal 

Figure  B.19.  Frequency  domain  of  recovered  signals,  layout  2. 
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Z  (deg)  |g  ,  Z  S,  (deg)  .g  ,  Z  S,  (deg) 


(a)  First  signal  (b)  Second  signal 

Figure  B.20.  Frequency  domain  of  recovered  signals,  layout  3. 


(a)  First  signal  (b)  Second  signal 

Figure  B.21.  Frequency  domain  of  recovered  signals,  layout  4. 


(a)  First  signal  (b)  Second  signal 

Figure  B.22.  Frequency  domain  of  recovered  signals,  layout  5. 
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(a)  Time-domain  signals  (b)  Time-domain  error 

Figure  B.23.  Time-domain  signals  and  error,  layout  1. 


(a)  Time-domain  signals  (b)  Time-domain  error 

Figure  B.24.  Time-domain  signals  and  error,  layout  2. 
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(a)  Time-domain  signals  (b)  Time-domain  error 

Figure  B.25.  Time-domain  signals  and  error,  layout  3. 


Time  (ms) 


(a)  Time-domain  signals  (b)  Time-domain  error 

Figure  B.26.  Time-domain  signals  and  error,  layout  4. 


(a)  Time-domain  signals  (b)  Time-domain  error 

Figure  B.27.  Time-domain  signals  and  error,  layout  5. 
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APPENDIX  C: 
Code  Listings 


C.l  delay  s_from_lay  do  wn.m 


function  [  tau  ]  =  delays_f rom_laydown {  L  ) 

%  DELAYS_FROM_LAYDOWN  Compute  the  delay  matrix  TAU  from  laydown  L 
%  Compute  signal  delays  from  a  set  of  transmitters  to  a  set  of 
%  receivers,  with  geometric  laydown  as  specified  in  L. 

Re  =  6371e3;  %  m 

num.recv  =  length (L . r . x) ; 
num.xmit  =  length (L . t . x) ; 

%  loop  to  compute  all  elements 
for  m  =  l:num.recv 

for  n  =  1: num.xmit 

theta  =  (L.r.x(m)  —  L.t.x(n))  /Re  *  180/pi; 

%  Law  of  cosines  to  get  from  theta  to  d(l,l) 
dx  =  sqrt  {  (L  .  r  .  alt  (m) +Re) '^2  +  Re''2  —  ... 

2* (L.r.alt (m) +Re) *Re*cosd (theta)  ) ; 
d(m,n)  =  dx  /  cos (atan2 (L . r . y (m) — L . t . y (n) , dx) ) ; 

end 

end 

%  Convert  from  distance  to  time,  using  an  approximation  of  the  speed  of 
%  light 

tau  =  d  ./  3e8; 

end 


C.2  block_deltas_from_delays.m 
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function  [  D,  cond_num  ]  =  block_cieltas_f rom_delays  {  tau,  omega  ) 
%BLOCK_DELTAS_FROM_DELAYS  Build  the  block  diagonal  matrix  of  delays 
%  [D,  cond_num]  =  BLOCK_DELTAS_FROM_DELAYS  returns  the  block  diagonal 

%  delay  matrix  using  the  propagation  delays  found  in  TAU  and  the 
%  frequencies  from  OMEGA. 

o, 

o 

%  The  condition  number  matrix  allows  checking  for  degenerate  cases 

[num.recv,  num.xmit]  =  size  (tau) ; 

N  =  length (omega)  ; 

D  =  exp (  — Ij  .*  reshape (tau, [num. xmit*num. recv  1])  *  omega); 

D  =  num2cell (reshape (D,  [num.recv  num.xmit  N] ) ,  [1  2 ] ) ; 

cond_num  =  squeeze (cellfun (@cond, D) )  ; 

D  =  blkdiag (D { : } ) ; 

end 


C.3  calc_signals_nb.m 


function  [  s,  r  ]  =  calc_signals_nb (  signal_spec,  t,  tau  ) 
%CALC_SIGNALS_NB  Calculate  sampled  narrowband  xmit  and  recv  signals 

o, 

0 

[num.recv,  num.xmit]  =  size  (tau) ; 

N  =  length (t ) ; 
dt  =  t  ( 2  )  ; 

%  Only  want  to  preserve  the  fundamental  freqs  and  magnitudes 
s  =  struct  (' mag '  ,  signal_spec  .mag,  'f,  signal_spec  .  f )  ; 

for  n  =  1: num. xmit 

s . func{n}  =  @ (t)  (s .mag (n) *cos (2*pi*s.f (n) *t)  )  ; 
s.time(n,  :)  =  s . f unc { n }  (t )  ; 

s.freq(n, :)  =  f f t ( s . func { n } (t ) ) *dt ; 

end 
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%  Received  signals  —  "r" 
r.time  =  zeros (num . recv, N) ; 
for  m  =  linum.recv 

for  n  =  l:num.xmit 

r .time (m, : )  =  r . time {m, : )  +  s . func {n}(t  —  tau{m,n)); 

end 

end 

r . sig_power  =  sum { r . time . ^2 , 2 ) /N; 
end 


C.4  recover_from_delays.m 


function  [  s_recovered  ]  =  recover_from_delays (  num,  D,  sig_freq,  dt,  . 
error_check;) 

%RECOVER_FROM_DELAYS  Use  signal  delay  to  recover  original  signal 

o, 

o 

%  Recover  the  original  signals  from  a  set  of  co— channel  signals 
%  recevied  at  a  separated  set  of  receivers,  given  the  known  delays 
%  between  the  transmitters  and  receivers. 

o, 

o 

%  s_recovered  =  recover_from_delays (  num,  D,  sig_freq,  dt, 

%  error_check  )  returns  the  original  signals  in  s_recovered.  D  is  a 

%  mN  X  nN  matrix  of  the  signal  delays  where  m  is  the  number  of 
%  receivers,  n  is  the  number  of  transmitters,  and  N  is  the  number  of 
%  discrete  frequencies.  sig_freq  is  the  matrix  of  frequency 
%  samples,  generally  the  result  of  f ft {time_domain_samples) .  dt  is 
%  the  sampling  interval  in  the  time  domain  and  error_check  is  a 

%  boolean  that  controls  an  extra  level  of  error  checking  in  the 

%  recovery  process,  useful  for  understanding  errors  at  discrete 
%  frequencies  due  to  poorly  conditioned  blocks  in  the  D  matrix. 

N  =  length (sig_freq) ; 

R  =  reshape (sig_freq,  [num. recv*N, 1] ) ; 

%  remove  the  rows  and  columns  that  correspond  to  omega (1)  =  0 
D  =  D (num. recv+1 : end, num. xmit  +  1 : end)  ; 
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R  =  R (num. recv+1 : end)  ; 


%  preallocate  S  for  speed  and  proper  dimension, 

%  necessary  to  account  for  the  missing  elements  per  above 
S  =  zeros {num. xmit* (N—1) ,  1); 

if  error_check 

%  cycle  through  the  blocks  to  ensure  independent  calculation 
%  there  are  N— 1  blocks,  after  the  removal  completed  above 
for  blk_num  =  1:N— 1 

index  =  2* (blk_num— 1)  +  ( 1 : num . xmit ) ; 

if  cond (D ( index, index) )  >  300 

fprintf('Bad  condition  number  at  indices:  [%d  %d]\n',  ... 

index+num . xmit )  %adjust  for  missing  block 
disp (D (index, index) ) 

end 

S (index)  =  D { index, index)  \  R (index); 

end 

if  {~all(D  *  S  —  R  <  eps)  |  -all  (D  \  R  —  S  <  eps)) 

warning (' Recovery  had  errors  greater  than  epsilon.’) 

end 

else 

S  =  D  \  R; 

end 

S  =  reshape  (S,  num. xmit,  []); 

%  insert  zeros  for  the  known  singularity  at  f=0  and  return  to  time  domair 
s_recovered . f req  =  [ zeros (num . xmit , 1 )  S]; 
s_recovered. time  =  if f t ( s_recovered . f req,  [],  2)/dt; 

s_recovered. sig_power  =  diag { s_recovered . t ime  *  s_recovered . t ime ' )  /  N; 

s_recovered . N  =  length (s_recovered. time) ; 

end 


C.5  calc_signals_wb.m 


function  [  s,  r  ]  =  calc_signals_wb (  signal_spec,  tau  ) 
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CALC_SIGNALS_WB  Calculate  sampled  wideband  transmit  and  receive 
signals 

[  s,  r  ]  =  CALC_SIGNALS_WB {SIGNAL_SPEC,  TAU)  returns  the 
transmitted  (S)  and  received  (R)  signals,  with  delays  specified  by 
TAU  incorporated  in  the  received  signals. 

SiGNAL_SPEC  is  a  structure  that  includes  SiGNAL_SPEC . TYPE,  which 
may  be  either  'noise'  or  'impulse'  for  a  AWGN  signal,  or  a  time 
domain  impulse,  which  has  a  constant  response  in  the  frequency 
domain.  The  impulse  is  useful  when  testing.  SiGNAL_SPEC . N  is  the 
number  of  points,  and  SIGNAL_SPEC . dt  is  the  inter— sample  interval, 
or  1/Fs. 

CALC_SIGNALS_WB {)  with  SiGNAL_SPEC . TYPE  =  'noise'  returns  a  signal 
with  power  spectral  density  equal  to  1. 

if  defined,  SiGNAL_SPEC . RNG_STATE  is  used  to  set  the  state  of  the 
RNG  prior  to  signal  generation.  This  ensures  consistent  signal 
contents  across  multiple  calls  to  CALC_SiGNALS_WB ( ) . 


s  =  signal_spec; 
dt  =  signal_spec . dt ; 

%  Use  a  seed  for  the  RNG,  if  requested 
%  constant 

if  isfield {signal_spec,  ' rng_state ' ) 
rng {signal_spec . rng_state) 

end 


keeps  signal  contents 


[num.recv,  num.xmit]  =  size  (tau) ; 

buffer_len  =  ceil {max (reshape {tau,  [num. recv^num.xmit  1]))  /  dt) ; 


for  n  =  1: num.xmit 

switch {lower {signal_spec . type) ) 
case  'noise' 

N  =  signal_spec . N; 

N_hat  =  2 . ^nextpow2 {N+buf fer_len) ; 
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s.time{n,:)  =  randn(l,N); 

s.freq{n,:)  =  fft {s . time (n, : ) ,  N_hat,  2)*dt; 

case  ' impulse ' 

N  =  signal_spec . N; 

N_hat  =  2 . ^nextpow2 {N+buf fer_len) ; 
s.freq{n,:)  =  ones(l,N_hat); 

s.time{n,:)  =  if ft {s . freq (n,  :  )  ,  N_hat,  2 ) / dt ; 
case  'ofdm' 

signal  =  ofdm_signal (s . ofdmParams { : } ) ; 

%  The  OFDM  signal  needs  to  be  interpolated  dt  time  domain 
N  =  ceil { signal . duration  /  dt) ;  %  new  number  of  samples 

N_hat  =  2 . ^nextpow2 {N+buf fer_len) ; 

t  =  ( 0  : N_hat— 1 ) -^dt ;  %  new  time  domain,  based  on  dt 

s.time(n,:)  =  interp_time (t ,  signal. sig,  signal . Tsamp) ; 
s.freq(n,:)  =  fft (s . time (n, : ) ,  N_hat,  2)*dt; 

end 

end 

df_hat  =  \  /  (N_hat*dt) ; 

fprint f ( ' s . time  energy:  %f\t%f\n',  diag(s.time  *  s.time')*dt) 
fprint f ( ' s . f req  energy:  %f\t%f\n',  diag(s.freq  *  s . f req ' ) *df_hat ) 

9-  9-- 
o  o 

%  create  a  new  omega  vector  for  the  longer  time  vector 

k  =  (0:N_hat— 1)  —  N_hat/2; 

omega  =  fftshift (2*pi*k  /  (N_hat*dt ) ) ; 

9-  9-- 
0  o 

%  Delayed  signals  (as  though  received  independently  at  the  receivers) 
for  m  =  l:num.recv 

for  n  =  l:num.xmit 

delayed . f req (m, n, : )  =  s.freq(n, :)  .*  exp (— 1 j *omega*tau (m, n) ) ; 

end 

end 

delayed. time  =  if ft (delayed. freq,  N_hat,  3)/dt; 
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for 


m 


1 : num . recv 


for  n  =  l:num.xmit 

del.t  =  squeeze (delayed . time (m, n, :)) ; 
del.f  =  squeeze (delayed . freq (m, n, :)) ; 
fprint f (' delayed . time ( %d, %d)  energy:  %g\n',  ... 

m, n, diag (del . t '  *  del.t)*dt) 
fprintf (' delayed. freq (%d, %d)  energy:  %g\n',  ... 

m, n, diag (del . f '  *  del . f ) *df_hat ) 

end 

end 

o,  g_ 

0  o 

for  m  =  l:num.recv 

for  n  =  l:num.xmit 

corr_plot { s .time (n, : ) ,  ... 

permute ( squeeze (delayed . time (m, n,  :))  ,  [2  1])); 

end 

end 

o,  g_ 

0  o 

%  Received  signals  —  "r" 
r.time  =  zeros (num . recv, N_hat ) ; 
for  m  =  l:num.recv 

for  n  =  l:num.xmit 

r . time (m,  : )  =  r . time (m,  : )  +  ... 

permute (squeeze (delayed. time (m, n,  :))  ,  [2  1]); 

end 

end 

% r.time  =  r . t ime ( : , 1 : N) ; 

r.freq  =  fft (r.time,  N_hat,  2)*dt; 

fprintf (' r .time  energy:  %f\t%f\n',  diag(r.time  *  r.time')*dt) 
fprintf (' r . freq  energy:  %f\t%f\n',  diag(r.freq  *  r . f req ' ) *df_hat ) 


9-  9- 
o  o 

%  Signal  powers 

s.sig_power  =  diag(s.time  *  s.time')  /N_hat; 
r.sig_power  =  diag (r.time  *  r.time')  /N; 
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end 


C.6  interp_time.m 


function  output  =  interp_time {  t,  x,  Tsamp  ) 

%  INTERP_TIME  Interpolate  using  the  ideal  low  pass  filter. 

%  INTERP_TIME (  tO,  x,  Tsamp  )  returns  an  interpolated  estimate  of 
%  x(tO)  where  tO  is  the  time  of  the  desired  point  for  interpolation, 
%  X  contains  samples  of  the  function  x{)  such  that  x[n]  =  x{n*Tsamp) 
%  and  Tsamp  is  the  sampling  interval. 

o, 

o 

%  INTERP_TIME (  t,  x,  Tsamp  )  returns  a  column  vector  of  interpolated 
%  values,  where  interpolation  is  performed  at  times  specified  in  t. 

o 

0 

%  Authored  by  Ned  Patterson 

n  =  0 : length (x) —1; 
n  =  ones { length (t ), 1 )  *  n; 

t  =  reshape (t,  [],  1);  %  Ensure  that  t  is  a  column  vector 

t  =  t  *  ones(l,  length (x)); 

%  Assume  we  have  sampled  at  the  nyquist  rate 
B  =  1/ {2*Tsamp) ; 

h  =  sine {  t/Tsamp  —  n  ); 

X  =  reshape (x,  [],  1);  %  Ensure  that  x  is  a  column 

%  Compute  the  sum  of  the  product  using  the  inner  product 
output  =  h  *  x; 
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